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Abstract 

In shallow-water waveguides a propagating field can be decomposed over three kinds of modes: 
the propagating modes, the radiating modes and the evanescent modes. In this paper we consider 
a semi-infinite waveguide with a free surface and an uneven bottom topography, and we analyze 
the effects produced on the coupling mechanism between these three kinds of modes. Using an 
asymptotic analysis based on a separation of scales technique we derive the asymptotic form of the 
distribution of the mode amplitudes, and the coupled power equation for the propagating modes 
describing the energy propagation. We show that the randomly perturbed geometry of the waveguide 
induces exponential decaying, with respect to the size of the perturbed boundaries, of the total 
energy carried by the propagating modes and we give an expression of the decay rate. 
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Introduction 

Acoustic wave propagation in shallow-water waveguides has been studied for a long time because of its 
numerous domains of applications. One of the most important application is submarine detection with 
active or passive sonars, but it can also be used in underwater communication, mines or archaeological 
artifacts detection, and to study the ocean's structure or ocean biology. Shallow-water environments 
are used to model acoustic propagation media such as a continental shelves, but these environments are 
very complex media because of indices of refraction with spatial and time dependences. However, the 
sound speed in water, which is about 1500m/s, is sufficiently large with respect to the motions of water 
masses that we can consider this medium as being time independent. Moreover, the presence of small 
spatial inhomogeneities in the water, heave or the ocean bottom roughness can induce significant effects 
over large propagation distances. 

The effects of random boundaries has been studied in many physical setups. For instance, in fluid 
mechanics, fluid flows in a medium with random boundaries has been studied in [5], for shallow- water 
wave propagation with a free surface or a random depth has been studied in [51 [13] , and also in wave 
propagation in underwater acoustics with a perturbed sea surface [71 1231 1211 127] . In [23] the author 
proposes new deterministic boundary conditions of Robin's type at the surface of the ocean using an 
expansion of the acoustic field at the water-air interface. These new boundary conditions are used in 
[21] to study wave propagation in a shallow-water waveguide with a free random surface. In [71 127] the 
authors use a conformal transformation of the perturbed domain to study the propagation problem. 
However, in both cases, the authors observe an attenuation and a phase modulation of the transmitted 
signal caused by the perturbations of the propagation medium. 

Mathematical studies regarding acoustic wave propagation in randomly perturbed waveguides have 
been carried out in many papers [TTl [T21 [131 113 EH ISl 120] , but only with the rigid-lid assumption at the 
waveguide boundaries, and with random perturbations inside the waveguide through variations of the 
index of refraction of the propagation medium. In this paper, we consider a two-dimensional shallow- 
water waveguide model with a Pekeris profile, and with a free surface and an uneven bottom topography 
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(see Figure[T]). Waveguide model with a free surface and an uneven bottom topography is a more common 
scenario in practice than the ones studied in |15lll7lED] under the rigid-Hd assumption. In our context, 
even if the random perturbations are small we study wave propagation over large propagation distances, 
so that the truncated expansion of the random boundary conditions considered in [23 to obtain Robin's 
type conditions cannot be used. In this paper, we consider a conformal transformation approach of 
the shallow-water waveguide model with a free surface. In fact, the conformal transformation is a 
convenient tool allowing to use the modal decomposition of the unperturbed shallow-water waveguide 
(see Figure [l] (a)). In this model a propagating field can be decomposed over three kinds of modes: 
the propagating modes which propagate over long distances, the evanescent modes which decrease 
exponentially with the propagation distance, and the radiating modes representing modes which can 
penetrate under the ocean bottom. Using an asymptotic analysis based on a separation of scale technique 
precisely described in Section [T] we show in this paper that the random perturbations of the waveguide 
geometry induce a mode coupling between the three kinds of modes. Theorem |4.1[ Theorem |4.2[ and 
Theorem |4.3| describe the coupling mechanisms between the propagating and the radiating modes in 
term of diffusion models but taking into account the effects produced by the evanescent modes. The 
coupling mechanism with the evanescent modes induces a mode-dependent and frequency-dependent 
phase modulation on the propagating modes, the coupling mechanism with the radiating modes, in 
addition to a mode-dependent and frequency-dependent phase modulation, induces a mode-dependent 
and frequency-dependent attenuation on the propagating modes. More particularly, we show in Theorem 



5.1 that the propagating modes lose energy in form of radiation into the ocean bottom and their total 



energy decays exponentially with the propagation distance. These phenomenon are summarized in 
Theorem 1 1 . 1 1 and constitute the main results of this paper. Let us remark that all the results obtained 
in this paper are also valid for electromagnetic wave propagation in dielectric waveguides and optical 
fibers ^22.Mi- 

Numerical simulation in waveguides with randomly perturbed boundaries is a very hard question, 
and can lead to very complicated methods to manage the medium geometry. To simplify this issue, 
several works have proposed simpler propagation models ^ i23i i24i |27, which can be easily managed 
by numerical methods. As a result, this paper proposes an asymptotic diffusion model in shallow- 
water waveguides with unperturbed geometry (see Figure [l] (a)), describing the propagation of the 
propagating modes which can be easily managed by numerical methods, and taking into account the 
random perturbations of the initial propagation medium (see Figure [l] (6)). 

The organization of this paper is as follows. In Section [l] we present the waveguide model, we state 
the main results of this paper, and we give a summary of the tools used to obtain these results. The 
remaining of the paper consists in introducing more precisely the tools, techniques, and some simplifying 
assumptions used to prove the main results. In Section |3] we study the mode coupling mechanism when 
the three kinds of modes are taken into account, and in the same spirit of [11, Chapter 20] and pTl [T^ . 
we derive the coupled mode equations. In Section |4j under the forward scattering approximation, we 
study the asymptotic form of the joint distribution of the propagating and radiating mode amplitudes. 
We apply this result in Section |5] to derive the coupled power equations for the propagating modes. In 
this section, we study the influence of the radiating and evanescent modes on the propagating mean 
mode powers. We show that the total energy carried by the propagating modes decays exponentially 
with the size of the randomly perturbed section of the waveguide and we give an expression of the decay 
rate. Finally, in Section|6]we first present the modal decomposition associated to unperturbed waveguide 
model (see Figure [T| (a)) and studied in detail in ^35]. Then, we prove Proposition 2.2 Proposition 2.3 
Proposition |3.2| and Theorem |4.1[ 



1 Waveguide Model and Main Result 

Throughout this paper, we consider a two-dimensional linear acoustic wave model In this context our 
propagation model is a planar waveguide. In shallow- water environment, modeling continental shelves 
for instance, sound propagation can be described by Euler's equations [HI Appendix A], so that the 
linearized conservation equations of mass and momentum are given by |11L Chapter 2] 

^ -dtp + V.u = 0, 



Ke{z,x) ■ ' (1) 

p^{z,x)dtU + \/p^ F, 
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Figure 1: Illustration of two semi-infinite waveguides modeling shallow- water environments with a 
Pekeris profile. In (a) an unperturbed waveguide with the rigid-lid assumption is considered. In (b) a 
perturbed waveguide with a free surface and an uneven bottom topography is considered. In the two 
cases the waveguide are homogeneous in the interior with propagation speed profile 0(2;, x) equal to ci 
in the ocean section of the waveguide, and cq in its bottom. Here, e is small parameter characterizing 
the amplitude of the random fluctuation of the free surface and the bottom topography (see Equation 
0). 

where p is the acoustic pressure, u is the acoustic velocity, pe is the density of the medium, and K^^ 
is the bulk modulus. The coordinate z represents the propagation axis along the waveguide, and the 
coordinate x represents the transverse section of the waveguide (see Figure [l] or Figure [2]). The forcing 
term F(i, z, x) is given by 

¥{t,z,x)^^{t,x)5{z^Ls)e,, (2) 



where is the unit vector pointing in the z-direction. The source profile '^{t,x) is defined by (21), so 
that "^(t^x) ~ f(t)S(x — xq). Consequently, the forcing term 

F(i, z, x) ~ fit)5ix - xa)S{z - Ls)e, (3) 

models a source which is closed to a point-like source located at (Ls^xq), and emitting a signal f(t) in 
the z-direction (see Figure [2] (6)). Let d > be the average ocean depth, we assume that the medium 
parameters are given by: 

K(zx) = l^' ,e(-oo+oo) 
i^,[z,x) xe [rf+V;^(z),+oo), 00,+ooj, 

( ^ / Pi if (y/(z),d + v;^(z)), ^( ^ . 

Pe[Z,X)^<. .„ ^ J , T/e^ ^ , \ Z e (-(X),+00), 

{ Po II X e [d+V^^(z),+(X)), 

where and Vjf model respectively the free surface and the bottom topography (see Figure [2] (a)). 
In our context, the free surface can model heave produced by the speed of the wind, and the uneven 
bottom topography can model a sandy bottom for instance with variations produced by water currents 

m- 

From the continuity of the pressure and the vertical component of the particle velocity across the 
transition from the ocean to the bottom, we have 



lim p(i, X, z) — lim p(i, x, z) 



and 



lim dTP(t,x,z)= lim dTP(t,x,z), 
so that we obtain from equations ([T]) the wave equation for the pressure field 
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Figure 2: Illustration of the shallow- water waveguide model. In both figures d represents the mean 
ocean depth, and the random fluctuations are given by the graphs of at the free surface, and d+V^ at 
the bottom of the waveguide model. In (a) we illustrate the medium parameters of the waveguide model 
and its random fluctuations at the free surface and at the bottom. In this flgure, the bulk modulus 
Kg{z,x) (resp. pe{z,x)) is equal to Ki (resp. pi) in the ocean section of the waveguide and Kq (resp. 
Po) in its bottom. In (b) we illustrate the source term F{t, z, x), which approximately models a point-like 
source at {Ls, xq) and emitting a wave f{t) is the z-direction. L/e characterizes the size of the section 
in which the random fluctuations are included (see (|7])). The arrow with a and going to the left at 
z = L/e indicates that no wave is coming from the right at the end of the random section. 



with velocity field 

c.M = VK.M/p.M^[l ;[ -e(-oo,+oox (6) 

and A — + (see Figure [l] (6) for an illustration). 

Moreover, from the continuity of the pressure field at the free surface x — V^{z), which is tanta- 
mount to neglect the surface tension, the wave equation (|5| is complemented by the following boundary 
conditions 

p{t, Z, V7(z)) = Psurface V(t, z) G [0, +Co) X (-CX), -|-Oo), 

where Psurface is the atmospheric pressure. However, one can assume without loss of generality that 

Psurface — 0; 

and then consider Dirichlet boundary conditions at the free surface, which corresponds to a pressure- 
release condition. 

In this paper we consider the Pekeris waveguide model ^ . This kind of model has been extensively 
studied for half a century , and has been widely used to model an ocean with a constant propagation 
speed profile (see Figure [T]) . Such conditions can be found during the winter in Earth's mid latitudes 
and in water shallower than about 30 meters [53] . The Pekeris profile leads us to simplified algebra but 
it underestimates the complexity of the medium. However, the analysis carried out in this paper can 
be extended to more general propagation speed profiles. This model could also be used to study the 
propagation of electromagnetic waves in a dielectric slab or an optical fiber with randomly perturbed 
boundaries. In fact, in the electromagnetic setting the bulk modulus and the density of the medium 
become the permeability and conductivity [551 1221 EH EH] ■ 

In the definition of the medium parameters Q and ([6]), and model respectively small fluctu- 
ations of the free surface and the bottom topography of the waveguide, and are given by 

V:iz) = ViVs{z)Uez), and V^iz) ^ V~^Vbiz)fbiez), with e « 1. (7) 

We assume that Vg and Vt are two independent mean zero stationary bounded stochastic processes, 
and fs and fb are two smooth functions with support included in (er],L) (with ry > 0), so that fsi^z) 
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and fb{ez) represent the localization of the random fluctuations of the free surface and the bottom 
topography of the waveguide over the section z G (rj, L/e). 

The scaling of the localization functions as been chosen according to the size of the random fluc- 
tuations {y/e ^ 1), so that we have to wait for long propagation distances (of order 1/e) to observe 
significant cumulative effects on the pressure wave (see (13) and Theorem For instance, wind 



speeds at about 5m/s induce a standard deviation of the surface roughness at about 0.1m, and the 
same order of magnitude can be considered for the roughness of sandy bottoms (see [HJ Chapter 2] and 
^^). Moreover, the correlation lengths Ic of the surface and bottom roughness in ([t]) are considered 
to be about 10m, that is large compared to the roughness magnitude 21, Chapter 2]. In fact, in most 
ocean-surface conditions the ratio between the standard deviation and the correlation length is small, 
and regarding the bottom roughness the correlation length is larger than its standard deviation (see 
[23]). Let us note that the two random processes Vs and H in ([t]) are considered to be independent. 
In our context, the reason is that the surface and bottom standard deviations are small compared to 
the ocean depth (d ^ 30m), so that the bottom topography does not produce any ripples on the free 
surface. Finally, for technical reasons we need to introduce regularity conditions for the free surface and 
bottom topography. First, let us note that we have assumed that the free surface is given by the graph 
of the process V^, so that the surface waves caused by the wind are not breaking. In fact, breaking 
waves cannot be described as a graph of a function. However, the derivative of a graph describing a 
wave which begins to break is not defined at the wave breaking point. In this paper, we assume that 
the trajectories of are of class with bounded derivatives, so that our free surface model describes 
a smooth ocean surface with small fiuctuations. We refer to [21, Chapter 2, Section 2.4] for an example 
of an autocorrelation functions allowing such a regularity for the free surface V^. Moreover, we also 
assumed that the trajectories of are of class with a bounded derivative, so that the bottom topog- 
raphy Vff models a smooth ocean bottom with small reliefs. We refer to 21, Chapter 2, Section 2.9] for 
an example of an autocorrelation functions allowing such a regularity for the bottom topography V^, 
but let us note that this assumption is physically relevant according to experimental measurements |2H 
Chapter 5]. 

Regarding the mixing properties of the free surface and the bottom topography V^, we assume 
that Vs and Vb in ([t]) are 0-mixing processes [55^ , that is considering 

J^z = J^o.z = criVsiu),Vb{u), 0<u<z) and J^z,+oo = criVs{u),Vb{u), z<u), (8) 

the (T-algebras generated by (Fs(w)j H(w))o<ti<2 and {Vs{u),Vb{u))z<u-, we have 

sup |P(A|B) -P(yl)| < 0(m), (9) 

z>0 

The (/)- mixing property ^ describes the decorrelating behavior of the random free surface Vs and bottom 
topography Vt, where e i^(M) n i^/^(M) is a nonnegative function characterizing the decorrelation 
speed. In other words the mixing condition ([9]), describes the decreasing dependence between and 
!Fz+z,+oo as Z — > +00, so that for Z large enough Tz and J-z+z,+oo are almost independent. The mixing 
condition ([9| is a technical requirement allowing us to prove the mains result of this paper (Theorem 

In this paper we propose an asymptotic analysis based on a separation of scale technique. The 
important scales of the problem considered in this paper are: The typical wavelength A of the pulse f{t) 
m (|3|, the propagation distance = L/ e, the correlation length Ic, and the standard deviation a^^ — ^fk 
of the surface and bottom roughness. These parameters are chosen so that we are in the well-known 
weakly heterogeneous regime [TT] : 

\-H^ - 1, \-^L^ > 1, and (Te < 1. (10) 



In ( 10 1, the first relation means that the typical wavelength and the correlation length of the perturbation 
are of the same order, so that the acoustic pressure wave produced by the source can interact with all 
the scales of variation of the medium perturbations. The second one means that we consider a high- 
frequency regime, while the last one means that the amplitude of the perturbations are small. For 
example, the typical value of the propagation distance is about 10km, the typical wavelength in the 
ocean A and the correlation length 1^ are about 10m, and the amplitude of the fiuctuations are about 
10cm [211 US]. In fact, as already remarked we have to wait for long propagation distances to observe 
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Figure 3: Schematic representation of the change of coordinate from the unperturbed domain Dq onto 
the perturbed domain Dy^ ■ 

significant cumulative efi^ects on the acoustic pressure wave when the fluctuations of the free surface and 
the bottom topography are small. Moreover, we also have 

Le > d and X^^d^l. 

The first relation means that the the ocean depth (d 30m) is small compared to the propagation 
distance L^, so that our scaling takes place in a shallow- water regime. The second one means that 
the ocean depth and the typical wavelength are of the same order, so that we have a fixed number of 



propagating modes {N{uj) ~ X~^d) in (12 1. 

To study the asymptotic behavior of the acoustic pressure wave in the perturbed waveguide domain, 
we use a change of coordinates through a conformal transformation [35^ introduced in Section [2] from 

Do — {u + iv £ C, such that u > O} 

onto 

Dy^ — {z + ix £ C, such that x > V^{z)^, 

and allowing us to study the pressure wave p{t, z, x) in a waveguide domain with a flat surface (see 
Figure [3]). As a result, one can focus our attention on pQ{t,u,v) defined by 

Po{t,u,v) ^ p{t,z{u,v),x{u,v)), 

which represents the pressure wave p{t, z, x) solution of ([s]) in a waveguide domaine with a flat surface, 
and where z{u, v) + ix{u, v) = ^{u + iv) is the change of coordinates from Dq onto Dy^ ■ Thanks to this 
conformal transformation the acoustic pressure wave po satisfies now Dirichlet boundaries conditions at 
the fiat surface of the waveguide: 

Po(t,u,0) = 0, V(t, u) G [0,+oo) X (-cx),-|-cx)). 

Moreover, the original pressure wave p{t, z, x) can be recovered from po{t, u, v) by inverting the change 
of coordinates: 

p{t, z, x) = po(i, Re{^^^{z + ix)), /m(<i>"^(z + ix))). 

We show in Section [2] that this change of coordinates transforms the random fiuctuations of the free 
surface to random fiuctuations in the interior and at the bottom of the waveguide (see Figure |4|. 
However, to study the effects produced by these perturbations, it is more convenient to consider the 
propagating wave in the Fourier domain. 

Throughout this paper the propagation model is a linear model, so that the pressure wave PQ{t, u, v) 
can be expressed as the superposition of monochromatic waves by taking its Fourier transform. Here, 
the Fourier transform and the inverse Fourier transform, with respect to time, are defined by 



g{oj) - / 9{t)e^"'dt, g{t) = 



27r 



%uj)e-''^'du;. 



(11) 
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In the new system of coordinates Dq, and in the frequency domain one can decompose the monochromatic 
pressure field po{u!,u,v) as follows, 

Po{u,u,v)^^Pj{u,u)(j).j{u,v)+j p^{u},u)(j)^{uj,v)dj+ p^{uj,u)<j3j{uj,v)dj, (12) 



3 = 1 







, radiating modes evanescent modes 

propagatmg modes 

where the modal decomposition corresponds to the one of the unperturbed Pekeris propagation model 
recalled in Section 6.1 (see Figure [T] (a)). Here k{uj) = uj/cq is the wavenumber, cq is the propagation 



speed in the bottom of the waveguide model, and N{uj) is the number of propagating modes defined 



by (83l. The coefficients pj(w, u) (j G {l, . . . , iV(w) }) and pj{uj,u) (V7 G (— 00, fc^(w))), defined by 



(87 1 and (|88|), represent the amphtudes of each kind of modes (propagating, radiating, and evanescent 



modes). Moreover, (j)j{uj, u) {j G {l, . . . , N{oj)}) and (p-yico, u) (V7 G (—00, fc^(w))), defined by (81 1 and 
[84] are the eigenelements of the Pekeris operator defined by ( [77| . Thanks to this decomposition, one can 
study the effects of the random fluctuations of the waveguide on the pressure field pa{uj, u, v) through the 



coupled mode equations (31 1 introduced in Section 3.1 In decomposition (12 1, the propagating modes 



can propagate over long distances, while the amplitudes of the evanescent modes decrease exponentially 
with the propagation distance. Finally, the radiating modes represent modes which can penetrate under 
the ocean bottom. 

Before introducing the main results of this paper, we need to consider some simplifications. First, 
the change of coordinate we consider, thanks to Assumption [T] and Proposition |2.1| in Section |2] is 
an approximation of the conformal transformation introduced above which is more convenient for the 
asymptotic analysis without changing the overall results of this paper. Second, we recall that we 
are looking for the cumulative scattering effects over large propagation distances because of the small 
variations of the medium fluctuations. As a result, because of its dissipative behavior, the evanescent 



part of the mode decomposition (12) can be neglected after long propagation distances. However, 
their influences in the mode coupling mechanism can be incorporated in the mode coupling mechanism 
between the propagating modes and the radiating modes as described in Section [3?2l Nevertheless, to 



obtain respectively well defined coupled mode equations (31 1 and to incorporate the contribution of 
the evanescent modes in the coupling mechanism between the propagating modes and the radiating 
modes, we need to consider the three following assumptions. We consider Assumption [2] in Section 



3.1 which consists in neglecting the contributions of the high evanescent modes (7 G (— cx),— 1/^)) in 
the mode coupling mechanism, and the contribution of the low evanescent and radiating modes as well 
(7 G Here, we have introduced a new parameter ^ such that 

e < ^ < 1, 

which means that both of these parameters will go to in the asymptotic analysis carried out in this 
paper, but e first, then ^. We consider a radiation condition for the evanescent modes (7 G (—1/^, —^)) 
given by Assumption [3] in Section |3.2| and meaning that the energy carried by the radiating modes 
decay as the propagation distance becomes large. We also consider Assumption |4] in Section [33] which 
is the forward scattering approximation and consists in neglecting the backscattering effects during 
the propagation. Moreover, we consider Assumption [5] in Section |4.1[ which consists in neglecting the 
mode coupling mechanism between the radiating and the evanescent modes. Let us remark that in 
[Itij . in which the random perturbations take place in the interior of the waveguide, the mode coupling 
mechanism between the radiating and the evanescent modes do not play any significant role. Finally, 
we consider Assumption [6] which corresponds to a nondegeneracy condition on the modal wavenumbers 
of the propagating modes, defined in Section |6.1| that is we assume that they are distinct along with 
their sums and differences. 

Consequently, with these assumptions, the pressure field po can be expressed as follows, Vu G [0,-L], 



Po 



where the scaling u/e corresponds the size of the perturbed section in our waveguide model (see Figure 
[2] (&)), and a^'^(a;, u) represents the forward mode amplitudes unique solution of (56 1. In ( [T3| we 



have introduced the modal wave numbers (/3i(a;), . . . , /3Ar((^)(i^)) of the propagating modes and (-y/^, 7 G 
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(^,fc^(cL)))) for the radiating modes, which are defined in Section 6.1 The main result of this paper 



regarding the asymptotic behavior of the mode amplitudes is given by the following theorem. 

Theorem 1.1 Under Assumptions 7|6 and the mixing conditions (|9|, the forward mode amplitudes 
{a'^'^{uj, •))cg(o,i); solution of the differential equation (56 I, converges in distribution 

as e to a diffusion process a^(w, .) with infinitesimal generator given by (64). Then, the 

stochastic process {a^{(jj, ■)){e(o,i) converges in distribution on C"'^^'^-' x L^{0, k^{uj)) as ^ ^ to a diffu- 



sion process a (ijJ, .) with infinitesimal generator = £^^0(1*) defined by (170 



Moreover, under the assumptions of Theorem\5.1\ we have 



lim — In 



N{uj) 

E 



E[|a-(^,i)l' 



-Aoo (w), 



(14) 



with 



A^iu;) ^ inf (( - (r-(..) + r^^''(..)) + (A^ (^) + AT{^)))X, (15) 



where 



5f = {xe M"("\ X, > Vje {!,..., 7V(c.)} and \\X\\l^r.,^, = (^,^)k«(., = l}, 



and r'='^(a;), V^'^iuj), A^/iu}), and A^/icj) are defined by i\6^, iGT}, and Mj^ with ^ = 



First, let us remark that the infinitesimal generator defined by (70) depends only on the N(uj) 



discrete coordinates, so that the radiating part of the mode amplitudes aP{u!, .) remains constant, and 
only the propagating modes exhibit a diffusive behavior. Moreover, according to its definition, this 
infinitesimal generator can be split in three part describing respectively the the coupling mechanism 
between the Af(a;)-propagating modes (^i^), the coupling between the propagating modes with the 
radiating modes {C2), and the coupling between the propagating and the evanescent modes (£3 ). More 
precisely, the coupling mechanism between the propagating modes themselves can be described as a 
conservative diffusion mechanism (no loss of energy) . The coupling mechanism between the propagating 
and the evanescent modes induces a mode-dependent and a frequency-dependent phase modulation. 
Finally, the coupling mechanism between the propagating modes with the radiating modes, in addition 
to a mode-dependent and a frequency-dependent phase modulation, induces a mode-dependent and 
a frequency-dependent attenuation described by (14). This last relation means that the total energy 



carried by the propagating modes de cays exponentially fast with respect to the size of the perturbed 
section of the waveguide. In Theorem 1 1 . 1 1 and as discussed in Section^ the matrices T'^'^{ijj) and T'^'''{uj) 
represent the energy transfer between the propagating modes, while the coefficients A^'''(a;) and A'^^{u)) 
represent the energy dissipation caused the mode coupling between the propagating and radiating modes 
and induced by the fluctuations of surface and the bottom topography. Consequently, the dissipation 
rate (15) takes into account the energy transfer mechanisms in addition to the dissipative effects. These 



phenomena are in agreement with the results obtained in [71 1241 127j in which the authors have observed 
an attenuation and a phase modulation of transmitted signals in such perturbed environments. 

As a result, in the asymptotic e <C 1^ ^ 1 one can approximate formally the statistical properties of 
the monochromatic pressure field pq by 



L 



N{uj) 



with 



aO(a;, L) = M,{oj, De^^^i^^"^) Vj G {1, . . . , 7V(c.)}, 



where the argument "tpj^io, u) describes the mode-dependent and frequency-dependent phase modulation 
of the forward mode amplitude aj{uj, u), while the modulus Mj{uj, u) describes the mode-dependent and 
frequency-dependent attenuation (14), which means that the total energy of the transmitted monochro- 



matic wave decays exponentially fast as the size of the perturbed section of the waveguide increase. 

The remaining of this paper consists in introducing the tools used in this paper and the simplifying 
assumptions (Assumptions T]|6 1 in a complete way to prove Theorem [1.1 [ through Theorem - 
l42l Theorem and Theorem [O 



4.1 



Theorem 
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2 Conformal Transformation of the Waveguide 



The free surface of our waveguide model (see Figure [2| is not convenient for the forthcoming asymptotic 
analysis based on a spectral decomposition of the acoustic pressure wave. One strategy consists in a 
truncated expansion of the pressure wave at the free surface [23] with respect to the amplitude of the 
surface perturbations {^/e <C 1) to obtain Robin's type boundary conditions. However, this strategy 
cannot be applied in our context since the propagation distance is too large (of order e^^ 1). In this 
paper, we use another strategy which consists in a conformal transformation allowing us to transform 
our waveguide model with a free surface to a waveguide with an unperturbed surface (see Figure [3|. This 
strategy has already been used in several contexts [Tj |H1 HH 112 and will appear to be more convenient 
to flatten a random surface over large distance. 

In this section we introduce a conformal transformation, which is a convenient tool since it allows 
to use the modal decomposition introduced in Section 6.1 for an unperturbed waveguide (see Figure [l] 



(a)). We will see in Section [s] that this transformation by flattening the surface induces deformations 
in the structure of the waveguide, that is deformations of the bottom and also induces variations of the 
index of refraction of the waveguide (see Figure |4]) . Consequently, these perturbations are responsible 
to a mode coupling between the three kind of modes (propagating, radiating, and evanescent) which is 
studied in Section [XT] 

To introduce the conformal transformation we consider the waveguide domain with a free surface as 
a subdomain of the complex plan C: 



Dv ^ \ z + ix £<C : x>V' 



where is defined by ([7| (see Figure |3| , and the domain corresponding to the waveguide with a flat 
surface is then given by 

£»o = {■" + w G C : w > o|. 

We refer to |35j for the basic properties of conformal transformations. In order to provide an explicit 
expression of this transformation (from Dq onto Dy,) we need to consider boundary conditions at 
infinity. We have assumed that the random perturbations of the waveguide geometry are localized in 
the interval z € (0,L/e), then the difference between the conformal map and the identity map goes to 
as Z — u + iv goes to +oo. The following proposition, which is a consequence of [35, Theorem 14.8], 
proves the existence and gives explicit formula of a conformal transformation from Dq onto Dy^ . 

Proposition 2.1 There exists a bijective conformal map <i> from Dq onto Dy^ defined by 

<i>(u + iv) = z{u, v) + ix{u, v), 

with 

z{u,v)^u / du — — r- — and x{u,vj^v / du— rr, ^, (16) 

V(u, v) such that u + iv Cz Dq. 

Let us note that we can extend the conformal transformation <f> as a homeomorphism from the closure 



of Dq onto the closure of Dy_, thanks to |351 Theorem 14.19], so that the derivation of formulas (16 1 is 
a classical result about harmonic functions which can be found in [I1|3T] for instance. 



Consequently, using the conformal transform $ introduced in Proposition 2.1 one can consider the 
pressure wave in the new system of coordinates: 

Po{t, u, v) — p{t, z{u, w), x(m, v)), V(t, u, v) G [0, +oo) x (— oo, +oo) x (0, +oo), 

where p is the solution of Therefore, thanks to the conformal transformation, pQ{t,u,v) evolves 
in the unperturbed domain Dq and satisfies Dirichlet boundaries conditions at the fiat surface of the 
waveguide: 

PQ{t,u,0) = 0, V(t, u) e [0,+oo) X (-oo,+oo). 

Moreover, the original pressure wave p{t, z, x) can be recovered from pQ{t, u, v) by inverting the change 
of coordinates as follows, 

p{t, z, x) = pQ{t, Re{^^^{z + ix)),Im{^~^{z + ix))). 
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However, Formulas ( 16 ) involve the term z{u, 0) which is not convenient for the forthcoming analysis. 
In fact, at the boundary w = it is difficult to use the inductive formula z{u, 0) — u—U{z(., 0))(m), where 
U is the Hilbert transform defined in Proposition |2.2| in order to obtain estimates in the forthcoming 
asymptotic analysis. In the following proposition we give an approximation of z(u,0) and then an 
approximation of x{u,0) to motivate Assumption [l] introduced below. 

Proposition 2.2 We have for every p e (2, +oo) 

z{u, 0) = u - ll{V,'){u) + e^-'^/PAl{u), (17) 

and then 

x{u,0) = V:{z{u,0)) 



V:{u)~U{V:){u)-V:{u) + e'/^-'/PAl{u), 



(18) 



with supg IjAf ||ip( 



2\\LP{ 



< C almost surely, and where C is a deterministic constant. Moreover, 



in (17 1 and (18), U stands for the Hilbert transform defined by 

1 f f{u) 



U{f){u) 



where p. v. denotes the Cauchy principal value. 



-p.v. 



du. 



We detail how to obtain ( 17 1 in Section 6.2 and ( 18 1 comes from ( 17 1 in addition to a Taylor expansion of 
V7 defined by (0. The term e^/^-^/P irT^lSl comes from the product between dV^ /du and ^^pA{{u), 
and all the remaining terms in the Taylor expansion have a smaller order of magnitude. Let us remark 
that we prove this approximation only in L^(IR) with p < oo for the reason that the Hilbert transform 
is not a bounded operator in L°°(M), and we chose p > 2 so that Al{u) in (17 1 and A|(m) in (18 1 do 
not play any significant role in the forthcoming analysis. In particular we need to have I — l/p > 1/2 
{p > 2) in (|17]) and 3/2 - l/p > 1 in 

Consequently, for the sake of simplicity in the proof of Theorem |4.1| given in Section |6.5[ let us 
consider the following assumption. 



Assumption 1 

and 



z{u,0) = u-U(y,'){u), 



x{u,o) = v:{u)-uiv:)iu)-v:{u). 



(19) 
(20) 



Thi s assumption give us an approximation of the conformal transformation $ introduced in Proposition 



2.1 



^by neglecting the contributions of e^^^/PA\{u) and e'^/^~^/P^2(")- However, this assumption will 

not change the overall results of this paper because the contributions of the correctors e^~^/PA\{u) 
(resp. e^/^~^/Pyl|(u)) are too small with respect to y/e (resp. e) to produce any significant effects. In 
Assumption [T] we keep only the terms which produce significant effects. 

We can remark that ( 19 1 illustrates the fact that z{u,0) is closed to the identity map (since is 
defined by ([7|) but where the small corrector comes from the conformal transform. The same remark 



holds for (20l, at the boundary of the domain Dq, x{u,0) is closed to the original free surface but 
also with a small corrector coming from the conformal transformation of the domain. 

Finally, we describe in the following proposition the asymptotic behavior of the real and imaginary 
part of the conformal transformation given by (16 1 under Assumption [l] 



Proposition 2.3 \/p > 2, we have 



sup \z{u,v) — 7i| + \x(u,v) — v\ < Ce^^^ almost surely, 

u, v>0 



and 



sup 

U JQ 



+ 00 



E 



\z{u,v) — up + \x{u,v) — v\' 



Proposition |2.3| means, as expected, that the conformal transformation is closed to the identity map 
since the amplitude of the surface fluctuations are small. The proof of Proposition |2.3|is given in Section 



where C is a deterministic constant independent of e. 
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3 Mode Coupling in Random Waveguides 



In this paper, the stochastic effects produced on the wave propagation is studied using the spectral 



decomposition of the unperturbed Pekeris operator defined by (77 1 and introduced is Section 6.1 In 



this section we introduce the mode coupling mechanism induces by the random free surface and the 
uneven bottom topography of the waveguide model (see Figure [2]) . 

The propagation model ([s]) considered in this paper is a linear model, so that the pressure wave 
p{t, z, x) can be expressed as the superposition of monochromatic waves by taking its Fourier transform 
(111, and then satisfies the time-harmonic wave equation 



9fp(cj, z, x) + dlp{u, z, x) + 



cl{z,x) 



p{u}, z, a;) = V • F{u:, z, x), 



where is defined by ([6]), and F{u), z, x) is the Fourier transform of the source term F{t^ z, x) defined 
by (U). 

Let us first determine the mode decomposition of the pressure field entering the perturbed section of 



the waveguide (0, L/e). According to the spectral decomposition introduced in Section 6.1 and described 
in Section [T] the monochromatic pressure field can be expanded as follows in the unperturbed section 



p{uj, z, x) 



N(uj) 
3 = 1 



Py{uj,z)4>y{uj,x)d'y+ / py{uj, z)4>y{Lu,x)d'y, 



propagating modes 



radiating modes 



evanescent modes 



where the mode amplitudes are defined by (87) andJSSl), and where (j)j{ui^u) (j £ {l, • . ■ ,iV(w)}) and 
0^(a;,it) (V7 € {— 00, k'^{(jj))), defined by (81 1 and |84[ are the eigenelements of the Pekeris operator 



defined by ( 77 1 . In this paper we consider a source profile ^' in ([2| admitting in the Fourier domain the 
following mode decomposition. 



^ (l}j{uj,xo)(j)j{u},x) + / (t)j{u},xo)cl)'riui,x)d-y 
^ J(_s,_5)u(5,fc2(„)) 



(21) 



where xq e {0,d). The meaning of the bounds S and ^ is that the source does not excite the high 
evanescent modes nor the low radiating and evanescent modes. However, S can be arbitrarily large 
and ^ arbitrarily small, so that "^{t^x) ~ f{t)5(x — xq) models a source profile closed to a point source 
emitting a signal f{t). 

Let us note that ^ is introduced for technical reasons as explain in Section |3.1| Moreover, let 
us remark that the bound S has been introduced so that ^'(a;, .) belongs to H — i^(0,-|-oo). This 
requirement is used to determine the expression of the pressure field entering the perturbed section of 
the waveguide using a coefficient identification allowed by the spectral decomposition of the Pekeris 
operator (77 1 in H. The source term ^ implies the following jump conditions for the pressure field 
across the plane z = Ls [HI HZ] , 



p{uj,x,L+) -p{uJ,x,Lg) = 
dzP{^,x,Ll,) - dzp{uj,x,Lg) = 0, 

so that the pressure field z, x) for z e {Ls, 0] is given by 



(22) 



P(w, z, x) 



N{uj) 

E 

i 



n(w) -o I \ 



^1/4 



s^v^^(/)^(a;,x)d7 



o k|l/4 ' 



(23) 



l''l^</)-y(a;,x)d7, Vx > 0, 
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and where 



(24) 



c 



7,0 



Now, to study the evolution of the pressure field in the perturbed section (0,L/e) through a mode 
coupling mechanism, we consider the new system of coordinates introduced in Section [2] This transfor- 
mation allows us to flatten the free surface of the waveguide domain, and then allows us to decompose 
the pressure field with respect to the propagating, radiating, and evanescent modes. 
The pressure field p{uj, z,x) in the new system of coordinates is given by 



Po{uj,u,v) —p{ijj,z{u,v),x{u,v)), y{u,v) E (— 00,+cx)) X (0,-|-c»), 



(25) 



where z{u,v) and x{u,v) are defined by (16) under Assumption [T] As a result, thanks to Proposition 

m 

(26) 



2.3 the monochromatic pressure field po{uj, u, v) satisfies the following time-harmonic wave equation 
dlpoiuj, u, v) + d^poioj, u, v) + k'^{uj)Je{u, v)nl{u, v)po{ui, u, v) = 0, 



V(u, u) e [0,L/e] X [0,+oo), with Dirichlet boundary conditions 

Pq{uj,u,Q) — Q, Vu e (— 00, -Koo). 



Moreover, in (26 1, the Jacobian and the index of refraction are given by 

2 



J^{u, v) 

and V(w,f) e (— 00,-foo) x (0,H-oo) 
ne{u,v) — CQ/ce{z{u,v),x{u,v)) - 



dux{u,v)^ + (^dvx{u,v) 



(27) 



(28) 



ni = co/ci > 1 
1 



if 
if 



ve{0,x-\u,d+V,%z{u,v)))), , . 
ve [x-\u,d + V^%z{u,v))),+oo). ^ 



According to the time-harmonic wave equation ( 26 1 , let us remark that the conformal transformation 
<i> introduced in Section [2] transfers the surface perturbations to the interior of the waveguide through 
the Jacobian Jj, and to the bottom of the waveguide through the index of refraction n^^u, v) defined by 
(29 1 (see Figure |4|. Moreover, the conformal transformation induces new perturbations on the bottom 
topography in addition to the initial perturbations V/^. 



Thanks to the Dirichlet boundary conditions ( 27 1 , the effects of the perturbations caused by the 
Jacobian Jj and the index of refraction n^(u,v) in (26 1 can be studied using the modal decomposition 
introduced in Section 6.1 and corresponding to an unperturbed waveguide (see Figure [T] (a)). 

3.1 Coupled Mode Equations 

Let us begin this section with two remarks. First, J^{u, v)ri^{u, v)pq{uj, u, v) does not necessarily belong 
to H = L^{0, +00) since the Jacobian blows up as v goes to 0. To correct this problem we need to 
introduce the following assumption (Assumption [2]) on the modal decomposition (12 1 of the pressure 
field po{Ld, u, v) defined by (25 1. 



Assumption 2 



N{u) 

Po{uj,u,v) = ^ pj{uj,u)(pj{uj,v) + 



p^{u},u)(j)y{uj,v)dj, 



(30) 



where the modes amplitude are defined by (|87| and ( 88 1 . This assumption consists in neglecting the 
contributions of the high evanescent modes (7 € {—oo,—l/^)) in the mode coupling mechanism, and 
the contribution of the low evanescent and radiating modes as well (7 G (— CiC))' that it leads us to 



12 




Figure 4: Illustration of the conforirial transformation on the waveguide with a free surface and an 
uneven bottom topography. Figure (a) and Figure (b) represent waveguides with a free surface. In 
these figures the vertical dashed lines mean that the medium parameters are not perturbed. In (a) the 
bottom is not perturbed, that is Vjf = 0, and in (6) the bottom of the waveguide is perturbed. Figure 
(c) represents the waveguide after the conformal transformation. In this figure the perturbed vertical 
dashed Unes mean that the medium parameters are perturbed. Moreover, the bottom is perturbed 
by the conformal transformation. However, if the bottom was already perturbed before the conformal 
transformation (Figure (6)), the bottom of the waveguide in Figure (c) is perturbed by the conformal 
transformation in addition to the initial perturbations of the bottom. 
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simplified algebra in the proof of Theorem |4.1| In such a decomposition, the radiating and the evanescent 
part are separated by the small band {—^,0 with ^ ^ 1. The goal is to isolate the transition mode 
between the radiating and the evanescent part of the spectrum of the Pekeris operator ( 77 1 given by ( 78 1 . 
The cut-off wavenumber —1/^ in the previous decomposition allows to obtain rigorous derivation of the 
following coupled mode equations more easily, but do not change the overall results obtained in Section 
|4.1[ Consequently, according to Assumption [2] the coupling mechanism of the spectral components of 
the pressure field will be studied in the space 

C^(") xL2((-l/e,-e)U(e,fc2H)). 

Finally, we assume that e <C ^ <C 1 so that we have two distinct scales. In this paper, we consider first 
the asymptotic e goes to and in a second time the asymptotic ^ goes to 0. 

Consequently, the mode amplitudes in ( 30 1 satisfy in the following coupled modes equation for 
u e [0,L/e]: 



1=1 



+ k\uj) C'',{uj,u)py{<jj,u)di = Q 



,2 



(31) 



d 

Pj{uj,u) + J p^{uj,u) + k'^{uj) ^ C^i{uj,u)pi{uj,u) 



where 



du^ 

1=1 

+ k^{uj) C^y(a;,u)py(a;,M)d7' = 0, 

r+oo 

C^g{uj,u) ^ / {J^{u,v)n'^{u,v) - n'^{v))(t)r{u!,v)(f)s{u!,v)dv, (32) 
Jo 

for (r,s) e {{l,...,N{uj)} x (-oo, fc2(a;)))2, and 

n{v)^\ -i-V'^i>l |f -e(0,rf), (33) 
^ ^ [ 1 if u e [d, +oo). ^ ^ 

Under Assumption [2] the random coefficients C"^(w,m), and then the previous coupled mode equations 
in C^(") X L2((_i/^^_^) u (^,fc2(w))), are well defined since 3Ki^, > 0, such that sup„| J,(m, d) - 1| < 
Ki^t/v^, and 

M^,v) -K2V Vre{l,...,iV(o.)}x(-l/e,-e)U(e,fc2(c^)). (34) 

Moreover, sup„|n^(u, v) — n'^{v)\ is bounded and has a bounded support with respect to the transverse 
variable v. 

In the previous coupled mode equations (31 ) the coefficients C"^(a;, z) represent the coupling between 
the three kinds of modes (propagating, radiating and evanescent modes). Let us remark that 

Jc(u, v)n^{u, v) — Ti^{v) = n^{u, v){Je{u, w) — 1) + (n^iu, v) — n^{v)), 

so that the first part of the decomposition induces random perturbations in the core of the waveguide, 
while the second part which is the difference of the index of refractions induces random perturbations 
at the bottom of the waveguide (see Figure |4]) . 

Next, we introduce the amplitudes of the generalized right- and left-going modes a{uj, z) and 5(w, z), 
which are given by 

p,{u;,u) = ^^(a,(..,u)e'ft(-)- + 6,(c.,u)e-^ft(-)"), 



d 

du 



p,{u,u) = zy)3>)(aj(c.,u)e*ft(-)« -S,(L.,u)e-'^^(-)«), 



^p^(a;,M) = i7i/4(a^(w,u)e'v^" - 6^(cj,u)e-'V7«^^ 
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so that 



Vj e {1, . . . , iV(ti;)} and almost every 7 g k'^{uj)). From now on, let us note 



(35) 



From pTj ) and (351, we obtain the coupled mode equation in x H'^ x 1/^, — ^) for the mode 
amplitudes (a{uJ, u),b{uj, u),p{uj, u)) : 



■12/ \ ^i^) r^e ( \ 



where 



a^(a;, m) 



^ ^ ^I^(^^(^^,),.(A-y.).+S.(.,u)e-<^'+^)") 



2 7^ ^i/4yi/4 



Ar(_ 

a^'(ti;, u)e 



*fc2(,,) /■^^(-) q^,(a.,z.) ^^^^ „,„.(vr-V7)« 



+ 6-y'(w, u)e 



d7' 



1/5 7 



1/4 



p^/(w, u)(i7 e 



57(w,zi) = fc2(c.) V ^'V; (SKw,")e*^'" + fo/(c^,u)e"'^'" 



p^(a;,u) + 7p^(a;,M) + \/eg^(w,M) = 0, 



+ fc2(w) 



^^ilT^^ (a7'(^.^^)e'^" + V(c.,w)e-V7'«) ^7' 



+ fc2(a;) / C:^y(a;,w)py(w,it)d7'. 



(37) 



(38) 
(39) 



(40) 



Let us note that in absence of random perturbations, the mode amplitudes a(cj,M) and b{u},u) are 
constant. Moreover, to complement this system, we need boundaries conditions. First, let us recall that 



the pressure wave entering the random section at z = is given by (23). Second, we assumed that no 
wave is coming from the right at z — L/e (see Figure [2] (6)). Consequently, we can complement the 
coupled mode system ( 36 1 , ( 37 1 , and ( 38 1 with the boundary conditions 



a(w,0) = ao'^(w) and [ - ] = 6^^^) in "H^, 



L 



(41) 
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where ag^(a;) and b'^£^{uj) are defined by (35), (87 1, and (88), for respectively u = and u = L/e. In 



Q (w) represents the initial amplitudes of the right-going propagating and radiating modes at 



(Iflf, a, 

u — while 6^'^(a;) represents the initial amplitudes of the left-going propagating and radiating modes 
at u — L/e. The following proposition shows that ag^(a;) is closed to the mode amplitudes of the 



pressure field ( 23 ) coming from the right homogeneous part of the waveguide and entering the random 
section at z = 0. As a result, ag''(££) is an approximation of ag(a;) defined by (24) in the new system of 
coordinates introduced in Section 2 The convergence of toward in the following proposition 

translates the fact that in the new system of coordinates we have assumed that no wave is coming from 
the right at z = L/e (see Figure [2] (fe)). 

Proposition 3.1 We have 



limE 

e-»-0 



= 0. 



This proposition means that the pressure wave po('^j u — Q,v) (resp. u — L/e) is closed to the pressure 
wave p{uj,z — 0,x) defined by (23) (resp. z — L/e). Pro posi tion 3.1 and gives us an approximation of 
the coefficient ag^(a;) in H^, which is used in Theorem |4.l| and Theorem 4.2 to determine the initial 

Proposition |2 . 3| and the 



conditions of the diffusion processes. This proposition is a direct consequence o: 
Parseval's equality associated to the spectral decomposition introduced in Section |6.1| since we have 



-t-oo 



E 



\po{u),u^,v) -p{uj,u^,v))\^ + \duPo{uj,u^,v) - dup{i^,u^,v))Y 

+ 00 



dv 



< C 



E 



\z{u^,,v) — Ut.\ + \x{u^,,v) — v\^ 



where is either equal to or L/e. Let us remark that this inequality holds because is out of the 
randomly perturbed section of the waveguide (see Figure[2|, so that we can use the variable v € [0, -l-oo) 
of the unperturbed transverse section also for the original monochromatic pressure field p{u}, z, x). 



3.2 Influence of the Evanescent Modes on the Propagating and Radiating 
Modes 

In this section we study the influence of the evanescent modes on the coupling mechanism between 
the propagating and the radiating modes. The main goal is to obtain a more convenient coupled 
mode equation involving only the propagating and the radiating modes, but taking into account the 
contribution of the evanescent modes. However, to obtain the following proposition we introduce the 
following radiation condition for the evanescent modes meaning that the energy carried by the evanescent 
modes decay as the propagation distance becomes large. 



Assumption 3 



lim 

>-|-oo 



+ 00 



l/e 



p^f{uj, u)4>-y{u!, v)d^ 



dv = 0. 



Consequently, thanks to Assumption [3] and following the proof of [XT', Section 4.3], we can rewrite (36 1, 
( [37| , and (jstj in a closed form in x H^. 



Proposition 3.2 We have 
d 



^^a{uj,u) = H°'|(w,M)(a(tj,u)) + 'H."'/^{uj , u) (b{u} , u)) 

+ G°5(a;,w)(a(a;,u)) -I- G^_j(a;, u) (5(w, w)) ■ 
+ A%u), 



r; 

R 



(w,u) 



a, is / 



(42) 



— 6(a;, u) = H^°^(tj, u) (a(a;, u)) + H^J'j(a;, z) (b{uj, u)) 

+ G'^'^^iiv, u) {a{io, u)) + G'^'^iiv, u) (b{u;, u)) 
+ B'iu), 



(43) 
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Vm G [0,i/e], where 



sup \\A'{u)\\n^ + \\B'{u)\\n^ <C sup (\\a{io,u)\\n^ + \\b{uj,u)\\n^) (44) 

ue[0,L/e] ^ ^ t. ue[0,L/e] ^ ^ ^ ^ 

for every p > 2. 

Let us make two remarks. The first one concerns the fact that we take p > 2 in order to have 



3/2 — 1/p > 1 in (44), so that the correctors A"^ and B"^ are small enough to do not play any significant 



role in the forthcoming analysis. Second, the proof of this proposition given in Section [6 . 4| follows closely 
the one of [I7i Section 4.3]. The only difference comes from error estimates near the boundary v = 



and the use of several truncated expansions in (32 1. Moreover, to obtain these uniform estimates we 
use Assumption [2] and the behavior (34 1 of the eigenelements introduced in Section 6.1 as w — > 0. 

In (|42| and (|43]), H"'^(w, u), H"*^(w,u), H^^j(a;,u), H^J'^(a;,M), G'^'^^{uj , u) , G^^^j(a;,u), G^^j(a;,M) and 
G''J^^{u!,u) are operators from H'^ to itself defined by: 



yie 



(45) 



H°,c,7(w,w)(y) = H^|'^_^(a;,u)(y) 



^l/4yl/4 



yie 



(46) 



Helj(^.")(y) = H""j^^(w, «)(?;) 



yie 



/ — ; 



Hfj_^(c.,K)(2/) = H-^(..,«)(2/) 



^1/4^1/4 



,e-»(\/'/'+v^)«d7' 



= G^^4j(w,w)(2/) 



ifc4(a;) 



_Af(w) 



//1/4 



' y-y"dj" dudj' 



(47) 



Gr,,,(a.,^.)(y) = G^^,,,(..,^.)(y) = ^ / / ^^XTC^ e-VI^II^I-^V^^^ 



yydj" dudj' 



(48) 
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/\^'\\u\-if3j{uj)u 



v"l/4 



X 



y^iid'y" dudj' , 



C':^ ,{u;,u) 



dud'y' , 



where C^{uj,u) = if u ^ [0,L/e]. 

The operators H°j(a;,M) and H"^(a;, u) represent the coupUng between the propagating and the 
radiating modes with themselves, while the operators G°'^(w, u) and G"j(w,m) represent the coupUng 
between the evanescent modes with the propagating and the radiating modes. Moreover, R°'^^(a;,u), 

Rg'j ^ (w, u), R^'^^(a;, u), and R^'^ ^ (w, m) represent the excitation of the evanescent modes produced by 
the source term on the propagating and the radiating modes through the random perturbations. These 
terms are defined by 

Kt,M = ^^j^M ^ ^ r' ^ll^0^,(^,,„)e-^(«-^^)dye-ft(-)", (49) 

2 J_s ^I3j{uj) 



Kt'i^^^) = K'^.'M = ^ / ^ ^7^|\;:r' '^7'(^.^o)e-VI-r'l("-^^)dye-^V7^, (50) 



'^'^ 4 J_yJ_s J ^mW\ (51) 



X (/>-y/'(w, a;o)e 



/|7"|(«+"-Ls)jy'^y 



7 7 



X (/)^'/(a;, a;o)e 



-Vl7"l(«+"-is)rfy'^y 



|7'||«|-iy7tt^~ 



(52) 



Equations ( 42 ) and ( 43 1 describe the coupHng process between the propagating and the radiating 



mode including the influence of the evanescent modes and the evanescent part of the source term. 
However, it is difficult to get good a priori estimates on 



sup 

ze[0,L/e] 



(\\aiuj,u)\\n^ + \\b{w,u)\\ 



to be sure that these residuals A"^ and B"^ in Proposition |3.2| are small even after large propagation 
distance L/e. One way consists in introducing the following stopping "time" 

L' = inf(L>0, sup \\a{u:,z)\\n^ + \H^.z)\\u^>^\, 

where a is such that 3/2— a > 1, in order to limit the size of the random section and control 
the size of the residuals and S'. However, will see in Section [373] that under the forward scattering 
approximation, we have VL > 

limP(L'= < L) = 0, 

e-i-O 

which means that with high probability the residual terms are small. 
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3.3 Forward Scattering Approximation 

The forward scattering approximation is widely used in the hterature, and consists in assuming that 
the couphng between forward- and backward-propagating modes is neghgible compared to the couphng 
between the forward-propagating modes. In this case, the amphtude of the left-going modes b should 
converge to as e goes to since we have assumed that no wave is coming from the right (see Figure [2] 
(6)). A rigorous derivation of the forward scattering approximation is technically complex. The main 
technical problem of a rigorous derivation is that the mode amplitudes (a{uj, z), b{u!, z)) is not uniformly 
bounded on [0, -L/e]. To correct this issue we could think that it suffices to use in a first step the stopping 
times corresponding to the first exit times of closed balls with respect to the norm || as described in 
[551 Chapter 11]. Unfortunately there are two problems. The first problem is that it seems not possible to 
show a limit theorem on C ( [0 , L] , (H^ , 1 1 . 1 1 -^^ ) ) . In fact , if the convergence holds on C ( [0 , L] , (H^ j 1 1 ■ 1 1 ) ) 
the energy conservation property ( 54 ) should be also valid for the asymptotic diffusion process according 



to the Portmanteau's theorem p^. However, this conservation property contradict the energy dissipation 
of the limit process observed in Theorem |5.1[ As a result, it seems to be more appropriate to obtain a 



hmit theorem on C([0, i], "H^j^) where H'^ ^ stands for H'^ equipped with the weak topology. Moreover, 
the second problem is that these stopping times are lower semicontinuous with respect to the topology 
of C([0,L],(-H^,||.||^^^)), but not with respect to the topology of C([0, where stands for 

H'^ equipped with the weak topology. 

However, one can formally show using the diffusion approximation theorem proved in [TT that the 
coupling between right-going propagating modes and left-going propagating modes involves coefficients 
of the form 

-t-oo 

R-ji{uj, u) COS ((/3/(w) -I- I3j{uj))u)du, 

while the coupling between two right-going propagating modes or two left-going propagating modes 
involves coefficients of the form 



+ 00 



Rji{u),z) COS - Pj{uj))u)du 



V(j, I) £ {l,..., N{uj)] , and where Rji € L'^iR) is defined by ([63]). Consequently, in our context if we 



forget these technical problems, according to [HI US] the forward scattering approximation should be 
valid in the asymptotic e ^ under the following assumption, which means that the Fourier transform 
of its w-autocorrelation function possesses a cut-off wave number. 

Assumption 4 



+ 00 



R,i{ij,z) COS {{Pi{uj)+ 13, {uj))z)dz^O V(j,0 G {l,...,iVM}' 



where Rji G L^(M) are defined by (63 1 



As a result, under these assumption, there is no mode coupling between the right-going and left-going 
propagating modes, but there is still coupling between the right-going propagating modes which is 
described in Section |4] Let us remark that the continuous part (0, fc^(w)) of the spectrum corresponding 
to the radiating modes, plays no role in the previous assumption. The reason is that the radiating 
modes and the evanescent modes play no role in the coupling mechanism between the right- and left- 
going modes as it has been shown in [TTJ [TS] . 

Consequently, according to the forward scattering approximation we only consider in what follows 



the simplified version of the forward coupled mode equation (42 1 and (43 1 



a{uj, u) = H°°5(w, u) (a(w, u)) + G^,^^{uj, u) {a{uj, u)) + ^(w, u) + R;;f " (w, u) + A'{u) (53) 



dz 



in ■H^,and on [0,L/e], where H^'^, G°°^, R"'|^^ and R,'^ are defined by (|45]), (|46|, (|47|, (|48]), (|49|, 
(l50|, (ISTI), (1521), and with Vp > 2' 



,3/2-1/p 



sup \\A%u)\\n^ <C- 



sup \\a{uj,u)\\u-; 

«e[0,L/e] 
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The evolution equation 
introduced in Section |3.1 
proposition. 

Proposition 3.3 



lim P(L^ < L) 

€-»-0 



53] ) is complemented with initial condition a{u!,0) — aQ^{uj), where aQ^(u;) is 
Moreover, under the forward scattering approximation, we have the following 



0, where U = inf ( L > 



sup ||a(w,'u)|l«- > 
«e[o,L/e] ^ e 



The proof of Proposition 3.3 follows exactly the one of [TSl Proposition 4.1]. Consequently, under the 
forward scattering approximation, one can adapt the results obtained in |T7| Section 4.2] to the system 
( 53 1 in order to derive the local energy flux conservation for the propagating and the radiating modes, 
that is V 77 > and according to Proposition |3.1| 



lim P ( sup 

^«e[o,L/£] 



|ja(w,u)| 



> ?7 



0, 



(54) 



which means that the amplitude a(w, u) is asymptotically uniformly bounded on [0, i/e] in the limit e — >■ 
0. The derivation of (54 1 is based on Lemma 6.1 Proposition 3.1 and Proposition 3.3 Consequently, 



using the Gronwall's lemma. Proposition |3.1| and Proposition |3.3| we have the following proposition. 
Proposition 3.4 We have V77 > 



lim I 

e-S-O 



sup \\a{uj,u) - ai(w, u)||^c. > lA = 0, 

«e[0,L/e] « ^ 



where ai{uj , u) is the unique solution of the differential equation 
d 



du 



ai{uj,u) 



H°'^(a;, u)(ai(a;, u)) + G°'^(a;, u)(ai(ijj, u)) with ai(a;, 0) = ag(a;) 



(55) 



As a result, thanks to Proposition 3.4 and [3, Theorem 3.1], we can focus our attention in what follows 
on the process ai{uj,u) unique solution of (55 1, which is studied in detail in Section |4] This proposition 
implies that the information transmitted to the propagating and radiating modes on the evanescent part 
of the source is lost during the propagation over large distances. 



4 Coupled Mode Processes 



In this section, we study the asymptotic behavior, as e — t' first and then ^ — > 0, of the statistical 



properties of the mode coupling mechanism (55 1 in terms of a diffusion process. First of all, let us define 



the rescaled process according to the size of the random section [0,i/e] 



This process is the unique solution on [0, L] of the rescaled forward coupled mode equations 



du 



{uj, u) = {uj, {a^^H^, u)) + -^G:^^ (u, (S^'« (w, u)) 



(56) 



in 'Hf, with initial conditions a'^'^(a;,0) = ao(^) defined by (24 1, and where H""^ and G^'^ are defined 
by (451, (|46|), (47l, (48l. Let us recall that, thanks to Proposition 3.4 (56) describes the evolution 



on [0,L] of the forward propagating and radiating mode amplitudes under Assumption [T]-[3] and the 
forward scattering approximation (Assumption |4]) . Consequently, under all these assumptions, (56 1 
allows us to describe the stochastic effects produced on waves propagating according to However, 
to describe these stochastic effects in terms of diffusion processes we introduce in the next section two 
more assumptions (Assumption [s] and Assumption [6| in order to use the ^mixing properties ^ of the 
free surface and bottom topography in a convenient way. 
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4.1 Limit Theorem 



This section presents the main results of this paper. The two following results are based on a diffusion- 
approximation result for the solution of an ordinary differential equation with random coefficients in 
a Hilbert space. However, according to the formula of the conformal transformation ( |16[ ) and the 
coefficients ( 32 1 , we cannot directly apply the asymptotic results obtained in [TTJ [T71 [15] since the 
mode amplitudes a'^'^{uj,u) are not adapted to the (/)-mixing filtration In fact, in the conformal 
map formula (16 1 all the random trajectory of the free surface is involved in ([32| because of the 
convolution. This fact lead us to technical difficulties in the proof of Theorem 4.1 (Section |6.5[ ) to 
manage the random coefficients (32 1, and to use the mixing properties of the random perturbations (|9|. 



To avoid this technical problem we consider the following assumption. 



Assumption 5 



C^^y{iu,u) = 0, V(7,7') e (-l/e,-6u(e,fc2(c.)). 



This assumption means that we neglect the coupling mechanism between the radiating and evanescent 
modes, they do not interact with each others. Let us remark that in |17j . in which the random pertur- 
bations talce place in the interior of the waveguide, the mode coupling mechanism between the radiating 
and the evanescent modes do not play any significant role. Consequently, in what follows we consider 



the evolution equation (56 1, where ^ and ^ are defined by (45 1 



He,5,7(w,w)(y) = 



Af(w) 

V 



(57) 



G"'^j(a;,u)(y) 



1=1 



yi dudj', 



/|7'||«|-»ftM« 



(58) 



and G^°^ 
Theorem 



4.1 



= for 7 e fc^((jj)), where we recall that C'^{u!,u) = if u ^ [0,i/e]. Moreover, to prove 
we consider the classical [TT| [T^ [20] nondegeneracy condition. 



and Theorem 4.2 



, Pn{lj){'-^)) defined in Section 6.1 are distinct along 



Assumption 6 The modal wavenumbers {/3i{uj), . 
with their sums and differences. 

Let us note that this condition is not necessarily satisfied for nonplanar waveguide model. However, 
in this case, Assumption [6] can be still valid if the statistical properties of the random fluctuations are 
rotationally invariant pU] . 

In the proof of Theorem |4.1 the random operator G°^ in (56 1 can be treated following the ideas of 



[22] as it has been done in |171 119]. In fact, the technique developed in [22] does not require any mixing 
properties. However, the random operator H°'^ requires to introduce an approximation of the random 
coefficients (32 1 to turn this formula in a more convenient form, and then to use the (/)-mixing properties 



of the filtration 



in Lemma 6.1 Section 6.5.1 



The approximations of the random operators H^'^ and G""! are given in Appendix 



In what follows, to simplify the presentation we will first focus our attention to the random pertur- 
bations at the surface of the waveguide, that is we assume that = 0. Afterward, we will take into 
account also the randomly perturbed bottom topography of the waveguide. 

The proofs of the results in this section are based on a martingale approach using the perturbed- 
test-function method. In a first step we show the tightness of the processes, and in a second step we 
characterize all the subsequence limits by mean of a well-posed martingale problem in a Hilbert space. 



4.1.1 The Random Surface 



This section is devoted to the asymptotic analysis of the forward mode amplitude a'^'^{oj,u) unique 



solution of ( 56 1 , and where the random perturbations of the waveguide are only due to the random 
surface, that is V? = 0. 
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Theorem 4.1 In the case — 0, under Assumptions^^^ and the mixing conditions VL > the 
family {a'^'^{ijj, •))eg(o i); unique solution of (56l, converges in distribution on C([0, as e — > to 

a limit denoted by a^{uj, •). Here "H^^ stands for the Hilbert space "H^ equipped with the weak topology. 

This limit is the unique diffusion process on T-L^, starting from a^{uj) defined by (24 1, and associated to 
the infinitesimal generator 

Cf\u) = criu) + C-;l{u) + C-;liu), (59) 

where 



2 _ 



and 



(60) 



-^3:1 = E (^) (^^■^^. - 



Let us note that according to the infinitesimal generator £^''' defined by (59 1, the diffusion phenomena 
happen on the support of fs which represents the locahzation of the random perturbations of the surface 
of the waveguide ([T]). 

Here, we have considered the classical complex derivative with the following notation: If w = vi+iv2, 
then 9„ = I {dy^ — idy^) and = ^ [dy^ + idy^). The coupling coefficients are given by the following 

notations: V(j, /) G {l, . . . , iV(w)}^ and j ^ I 

^7^''^ = 0R( \r( ^ / RA^.z)cOB{{Pi{^)-P,{iv))z)dz^ 
ZPj(UJ)Pl(UJ) Jq 



1=1 



^y(^) = 9^r \r( ^ / fi.K^,^)sin((A(o.)-/3,(a.))z)dz, 



(61) 



1=1 



and V(.?,0 e {!,..., iV(c^)}' 



U4/, ,\ /■+00 



rk^(<^) i-4/ \ r+oo 

Af'«(c.)= / ^ , / i?,v(^,^)cos((v/y-/3,(a.)).)dzd7', 



2Vy/?,(c.) Jo 

1,4/ \ /•+00 ^ 

A^/'Hoj)^ ^ , / i?,y(^:^)sin((v/y-/3,(c.))z)d^dy 



(62) 



^^%)= I 'iP^ pG;^(..,z)cos(/3,(..)z)eVM.rf,rf^,^. + !^Gf 
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Here, we have 



G«(Z) =4K-l)202(rf)^2(^)_ 



16 



Rs{z + {w - w')d) , , , 

^ aw aw 

{w^ + d^){w'^ + d^) 

Rs{z + wv — w'v') 



(w2 + v'^){w'^ + V'^) 



dwdw'dvdv' 



0j('u)(/)^('u)(/)j('u')07(^'')« V 



8(ni-l) 



TT 

16(ni - 1) 

TT 

16 



(t)j{d)(t)^{d) J (t)j{v)<j)^{v) ^ 



+ wv — w'v') 

{W^ + w2)2(ti;'2 + ^;'2)2 

+ WV — w'd) 



dwdw' dvdv' 



(i)j{v)cj)^{v)cj)j{v')(i)^{v')v'^ 



[w^ + v^){w''^ + d^) 
Rs{z + wv — w'd) 

Rs{z + uiu — w'v') 



(^2 + t;2)2(u;'-^ _|_ yl-^-^ 



dwdw' dv 

dwdw'dv 
dwdw'dvdv' , 



and 



6rf3 



Rs{w - w') , , , 3c? 

^ dwdw 

{w^ + d^){w'^ + d^)^ 7r2 



+ 4(ni-l)0,(d)0'(d) 



2wZi(i?,)(0) 



2^2 



(u;2 + d2)(y,/2_^^2) 

ww' Rs{w — w') 

{w'^ + w2)2(y;2 ^ ^2)2 



(ii;2 + d2)(y,/2^^2) 

dwdw' 



dwdw' +U{Rs){0) 



dwdw' + -U{Rs){0) 

V 



1 



(t«2 + v^){w' + u2) 

, , , iv^ 



dwdw' 



(w2 + -i;2)(u;'2 + ^2) 



Rs{w-w') 



(w2 + i;2)2(y;/2 _|_ ^2) 



1,2^2 



-dwdw' 



In the previous expressions is defined by (33 1, Rg stands for the autocorrelation function of the 
random process Vs , and Rg is the autocorrelation function of the processes Vg and ^ Vs , which is j^Rs- 
Let us note that Rg is an even function, so that Rs{0) — 0. Therefore, we need to assume that Rg 
decrease to fast enough so that its Hilbert transform at 0, U{Rg){0), is well defined. 
Moreover, Rji{uj, z) and Rji{uj, z) can be split in 6 terms 

6 

n=l 
6 

n=l 

which represent the different coupling mechanisms caused by the perturbations of the waveguide. We 
have 



R 



(1) 



R 



(2) 



R 



(3) 



R 



(4) 



R 



(5) 



R 



(6) 



4d2(ni - 1)2 



cj,'id)q^Ud) 



Rs{z + w — w) 



-dwdw. 
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4>] (w, v)(j)i{uj, v)(j)j (w, v)(j)i{u}, v) 



(w2 + d'^)[w'2 + ^2) 

i?s(2 + w — li) 
(it;2 + v^)(w^ + 52) 

i?s(z + w — w) 



4>j{uj, v)4>i{uj, v)4>.j{uj, v)4>i{uj, v)v'^v^ 



(ni - 1) 



(j)j{uj,d)(f>i{uj,d) / <pj{uj,v)<pi{uj,v) 



md{ni ~ 1) 



(t>j{u,d)(t>i{uj,d) / u)0i(w,t;)v 



16 



v)4>i{uj, v)4>.j{uj, v)4>l{uJ, v)!)"^ 



(u)2 + v'^)'^{w'^ + u2)2 
i?s(2 + W — w) 

(^2 + d'^){w'^ + i;2) 
R8{z -\- w — w) 

(w2+d2)(li)2 + „2)2 
+ ^ 'J') 



dwdw dvdv, 

dwdw dvdv, 
dwdwdv, 

dwdwdv. 



(63) 



(l(j2 + w2)(tt;2 + {i2) 



zdwdwdvdv. 
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and 



5(2) 



r(3) 
^3l 



p(5) 
^3l 



R 



16 

8d(ni - 1) 



0j (cj, v)(l)j{uj, v)(l)i{uj, v)(t>i{uj, v) 



Rs{z + w — w) 



dwdwdvdv. 



(l)j{uj,d)(f)j{uj,d) / 0,(w,u)0i(a;,'i;) 



16d(ni - 1) 



(j}j{u},d)(j)j{u},d) / 



i?s(2 + w — tZi) 

(w^ + V^)^{'W^ + u2)2 

i?s(z + w — w) 

(w2 + d'^){w^ + u2) 

i?s(z + W — t«) 
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(w2 + d2)(tt,2 _|_ y2)2 
+ W — ?«) 



dwdwdvdv, 
dwdwdv, 

dwdwdv, 



{w^ + v^)(w'^ + v^y 



: dwdwdvdv. 



Here, R^-^"^ and i?'^-* represent the coupling mechanism between the propagating modes at the bottom 
of the waveguide, while R'^^\ R^'^\ R^^\ R^-^^ and Iv-^^ represent the coupling mechanism through 

the transverse section of the waveguides. Finally, R^'^\ R^^\ and i?*-^-* represent the coupling 

mechanism of the propagating modes between the bottom and the transverse section of the waveguide. 

Finally, let us remark that the convergence in Theorem 4.1 holds also on C([0,L], ("H^, ||-||«^)) only 
for the Af(a;)-discrete propagating mode amplitudes. 

4.1.2 The Random Surface and Random Bottom 

In this section we present the asymptotic analysis of the forward mode amplitudes a'^'^{u!,u), where we 
consider a randomly perturbed bottom topography in addition to the randomly perturbed free surface. 
This section presents no additional difficulties since the random perturbations of the surface and the 
bottom are assumed to be independent, and then can be treated directly using the method of pTl [TU] . 

Theorem 4.2 Lender ^ssMmpfionsj7||6j and the mixing conditions ([9|, VL > the family (a'^'^{uj, ')) ^^(^q ly 
unique solution of (56l, converges in distribution on C{[0, L\,'H'^ ^) as e ^ to a limit denoted by 
a^{ui, ■)■ Here H'^ ^ stands for the Hilbert space "H^ equipped with the weak topology. This limit is the 
unique diffusion process on Tif , starting from aQ{uj) defined by (24 1, and associated to the infinitesimal 



generator 



where >C^'*(u) is defined by (59), and 



with 



{ [u) = (u) + C^ ^ [u) + £3;^ (u) , 



(64) 
(65) 



3:1 = 1 



E r'i'H {T3TidT,d^ + T,Tid^dT,-T,TidT,dT,~T,Tid^d^ 



3:1 = 1 



+ i 



3 = 1 



and 



2,C 



N{uj) 



3=1 

Af(w) 



(66) 



i=i 
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Let us note that according to the expression of the infinitesimal generator and C^'^ defined 
respectively by ( 59 1 and ( 65 1 , the diffusion phenomena happen on the support of fs and /f, which 
represent the localization of the random fluctuation of the surface and the bottom of the waveguide Q . 

The coefficients r'='^'«(a;), V^^'^iuj), V^'^iuj), K^-^^^iuj), A''''''«(w), and K^^«(a;) are defined as follows: 
V(j,Oe {!,..., A^M}' and J 



Rb{z) COS {{Piiuj) - I3j{uj))z)dz, 



1=1 



2k^{u){ni~\fct>']{io,d)(t)'i{io,d) f+ 



(67) 



i?6(z)sin {{I3i{uj) - fij{uj))z)dz, 



1=1 



and V(j,0 e {!,..., N{uj)}\ 



2fc4(w)(ni - l)^(l)''j{uj,d)(j>f{uj,d) 



+00 



Rb{z)dz, 



(w) = I 



2fc4(w)(ni - l)202(^^rf)02^(^^rf) 

« 2fc4H(ni^l)202(^^rf)0 2^(^^rf) 

-1/? " 



+00 



-R^l-s^) COS (( Vt ~ Pj{u;))z)dzdj\ 



+00 



+00 



+ »-R6(0) 



8fc2(c.)(ni-l)202(^^rf)0/2(^^^) 

/3,(w) 



i?6(z) sin ((VY - /3j(w))z)dzd7', (68) 
i?6(z) cos (/3J(a;)z)e~^/l'^l^dzd7yJ■ 



The drift {C^''^, 'Cg'^, 'C2ij' ''-'Sf defined by (60l and (66l) and the initial value 0,^(10) defined 
by (24 1 of the diffusion process a^{uj,-) still depend on the parameter ^ introduced in Assumption [2] 
Before discussing the meaning of each terms in ( 59 1 and ( 65 ) , we obtain the following result regarding 
the asymptotic ^ H> 0. 

Theorem 4.3 VL > 0, the family (a^{uj, '))^^(^q converges in distribution on C{[0,L], {Hq , W-Ww^^)) as 
~> to a limit denoted by a,{uj, •). This limit is the unique diffusion process on T-Lq , starting from 

ao{uj) — lim ag(w), (69) 



where ao(w) is defined by (24), and associated to the infinitesimal generator 



C-{u) = C-^{u) + C^{u) + C-^{u), 



(70) 



where 



N{uj) 



\ J2 /.^(«)((Af (c.)+zA;'^(^))r,9T, + (Af (.;)-*Af (c.))T,% 
i=i 

+ /2(u) ((Af (c.) + ^Afiu;))T,^T, + (Af (..) - ^A^■''{u:))TJ^^ 



N{uj) 
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withyj e {!,..., N{oj)} 



A'r-^uj) = lim A^'^'^M, M = lim Af ^^(o;), k^^{u^) = frnn^/Ho^), (71) 



and 



Af{oj) = lim Af Af{Lj) = lim Af = lim ^^'^(w) 



5^0 



(72) 



(841, (851, and (86l 



Let us note that in Theorem 4.3 the hmits in (71 ) and (72) with respect to ^ are well defined thanks to 



Consequently, Theorem |4.1[ Theorem |4.2| and Theorem |4.3| describe the asymptotic behavior, as 
^ first and then ^ — > 0, of the statistical properties of the forward mode amplitudes a'^'^{uj, •), in 



terms of a diffusion process with infinitesimal generator ( 70 ) and starting from ( 69 1 . 

The infinitesimal generator is composed of three parts which represent different behaviors on 
the diffusion process. However, we can remark that this infinitesimal generator depends only on the 
Af(a;)-discrete coordinates, so that the radiating part of the forward mode amplitudes remains constant. 
This result has already been obtain in ^TTi [TO] in a different setup. The first operator £^ describes 
the mode coupling between the Af(a;)-propagating modes. This part is of the form of the infinitesimal 
generator obtained in [TTJ [TS], and for which the total energy is conserved. The second operator £2 
describes the coupling between the propagating modes with the radiating modes. This part implies a 
mode-dependent and frequency-dependent attenuation on the A^(a;)-propagating modes that we study 
in Section [5] and a mode-dependent and frequency-dependent phase modulation. The third operator 
£3 describes the coupling between the propagating and the evanescent modes, and implies a mode- 
dependent and frequency-dependent phase modulation. The purely imaginary part of the operator 
does not remove energy from the propagating modes but gives an effective dispersion. 

The main consequence of these limit theorems is the attenuation of the energy carried by the prop- 
agating modes. According to the formulas of A'^''' and A^'** defined in (62 1 and (68 1, this attenuation 



is induced by a mode coupling between the propagating and the radiating modes depending on the 
localization of the random perturbations. The mode coupling describes by A'^''' is relatively simple and 
is localized at the bottom of the waveguide, while the mode coupling mechanism caused by the random 
surface is more complex. We have seen that the random surface produces, through the conformal trans- 
formation (16 1, perturbations in the structure of the waveguide (see Figure HI), that is at the bottom 



of the waveguide and the core of the waveguide. Then, the random fiuctuations of the surface produce 
a mode coupling, described by A'^'*, between the propagating and the radiating modes at the bottom 
and also all along the transverse section of the waveguide. The dissipation phenomenon on the energy 
propagation is described more precisely in Section |5] 

In absence of random perturbations at the surface, just with random perturbations at the bottom 
of the waveguide, for lu large enough the diffusion effects on the low propagating modes (j ^ N{uj), 
where N{uj) is the number of propagating mode (83l) is very small compared to the one on the high 
propagating modes. More precisely, for all bounded smooth function (p we have 

V77e [0,l[,Vje{l,...,iV(c.)"}, lim {u)^{T,) = 0. 



This result comes form (81) and the fact that V77 e [0, 1[ we have 



lim sup \aj{uj)—jn\ — 0, 

j€{l,...,N{uj)v} 



where crj{u}) is defined by (82). In other words, for uj large enough the coupling mechanism of the 
low propagating modes with the other modes is very small, so that these modes almost do not see the 
random perturbations of the bottom of the waveguide. 



4.2 Mean Mode Amplitudes 

In this section we describe the effect of the random perturbations on the forward mean mode amplitudes. 
From Theorem 4.3 we get the following result. 
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Proposition 4.1 Vz e [0,L], Vj G {l, . . . , N{uj)} 



lim limE 



exp 



X exp 




f1{u)du 



fl{u)du 



where ao(w) is defined by (69) 



First, let us note that the mean amplitude of the radiating part remains constant on L^(0, fc^(ijj)), 
since the diffusion process only holds for the propagating modes (see Theorem 4.3 1. Second, Vj G 
{l, . . . , N{uj)], the coefficients {Tjj{uj) +A'j{uj) — Tjj{oj))/2 are nonnegative according to the Bochner's 
theorem [55] , 

The decay rate for the mean jth-propagating mode is given by 



E 



which depends on the effective coupling between the propagating modes, and the coupling between the 
propagating and the radiating modes. This decay describes the loss of coherence of the transmitted field 
caused by the random perturbations localized by fs and /„. However, according to [111 115], in which 
the study takes place in a waveguide with a bounded cross section so that the spectral decomposition 
has no continuous spectrum, the authors obtain the same result without the terms A'^''* and A'^'*'. We 
will see in the next section that these coefficients are very important in our context. In fact, in the 
case of waveguides with a bounded cross section the authors observe decoherence of the propagating 
mode amplitudes but with conservation of the total energy carried by these modes. However, as already 



observed in ^ITJ [12] and as we will see in the next section in Theorem 5.1 the total energy carried by 
the propagating modes is not conserved in semi-infinite waveguides. This loss of energy comes from 
the radiative loss (coupling with the radiating modes) of the propagating modes in the bottom of the 
waveguide (see Figure Isl. 



5 Coupled Power Equations 

This section is devoted to the analysis of the asymptotic mean mode powers of the forward propagating 
modes. More precisely, we study the asymptotic effects of the coupling between the propagating modes 
with the radiating modes to lighten the radiative loss, that is the loss of energy of the propagating 
modes through the radiative modes. Let us consider for j e {1, . . . , N(ll))} 



Tiiuj, z) — lim lim E 



(w,z)| 



= E 



|aj(u;, z) I 



be the asymptotic mean mode power of the jth propagating mode. In other words, 7}(w,z) is the 
expected power of the jth propagating mode at the propagation distance z. From Theorem |4.3| we 
obtain the following coupled mode power equations: Vj' € {!,..., N{uS)} 



d 
dz 



TjiLJ,z) 



^ {f^{z)A':'%u;) + f^{z)Af{u;))T,{u;,z) 

^ {f!{z)T';;;{uj) + f!{z)T:;^{io)) (r„(^, z) - z)) , 

n=l 



(73) 



;.o(w)| , and where ao{^) is defined by (69 1. These equations describe 



with initial conditions 7j(w, 0) = 

the transfer of energy between the propagating modes, where r'^'*(a;) and F^^o;) are the energy transport 
matrix defined by (61 1 and (|67|. Compared to the couple power equations obtained for waveguides with 



a bounded cross section [ITJ[TS], we also have the radiative coefficients A^''(w) and Aj''' (uj) defined by 
( 62 1 and ( 68 1 resulting from the coupling between the propagating modes with the radiating modes. 
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Figure 5: Illustration of the radiative loss in the shallow-water random waveguide model. Figure (a) 
illustrates the radiative loss of the propagating modes in the bottom of the waveguide caused by the mode 
coupling with the radiating modes. Figure (&) illustrates the energy conservation of the propagating 
modes. 



These radiative coefficients are responsible for the radiative loss of energy of the propagating modes 
in the ocean bottom (see Figure [s]) , that is the total energy carried by the propagating modes decays 
exponentiall y w ith the size of the randomly perturbed section. This loss of energy is precisely quantify 

'(w) + A''-'^{uj) — 0, it has been shown in [llj and ^ISj 



in Theorem 



5.1 



If there is no radiative loss A'^'^( 
Chapter 20] that the energy of the propagating modes is conserved and for large random section the 
asymptotic distribution of the energy becomes uniform over the propagating modes. 



Let us defined 



S 



{ 



X e 



X,>0 



Vj e {!,..., 7V(c^)} and 



with (X,y)^„,^, = Efir' X,Y, for {X,Y) e (M^("))2, and 



A^,{u) = AYico) + A^H = dtag{ArH + Al^\oj), A-^(c.)(c.) + A^''(a.)(c.)) 



(74) 



To simplify the presentation of the following theorem, let us assumed that the localization functions 
fs and /b of the random perturbations are given by the indicator function 1(q^)(z). 

Theorem 5.1 Let us assume that the energy transport matrix r'^'*(w) +r'^'''(a;) is irreducible and that 
the localization functions fs and ft of the random perturbations are given by the indicator function 
l(o,L)(z). Then, we have 



lim — In 

L— ^-t-oo Ij 



N{lu) 

E 



T,{u;,L) 



= -Aoo(l^) 



with 



A^iuj) = inf (( - (F^'^Xo.) + r^^^(c.)) + AS(c.))X,X) 



xes 



where the energy dissipation rate (75 I is positive as soon as one of the coefficients A 
positive. 



(75) 



iLo)+Afiuj) is 



The proof of Theorem 5.1 follows exactly the one of [T71 Theorem 6.1] and comes from [6 . This 
result means that the total energy carried by the expected powers of the propagating modes decays 
exponentially fast with the size of the random section, and the decay rate can be expressed in terms 
of a variational formula over a finite-dimensional space. Let us remark that the energy dissipation rate 



( |75| depends on energy transport matrix F'^''*(ci;) and r'^'''(a;) defined by ( |61| and ( |67| in addition to the 
dissipation coefficients A^{uj) defined by (74 1. Then, the energy dissipation effect takes into account the 
energy transfer mechanisms in addition to the cumulative dissipative effects caused by the fluctuations 



of the free surface and the bottom topography. Formula (75) of the dissipation rate Aoo(w) is not simple, 
but it has the following lower and upper bounds. 



j<£{l,...,N(u,)} 



A^'^(t^)+Af (w) < A^iio) < A(w) 



1 



N{uj) 

E 



A^^'{co)+A^'\u;). 



(76) 
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In what follows we discuss the dissipation rate ( |75| according to the energy transfer mechanisms de- 
scribed by ^'^{uj) and r^'^(ijj), produced respectively by the fluctuations of surface and the bottom 
topography, and give an explicit expression of Aqo (w) • We will see below that the stronger the energy 
transfer is the larger the dissipation rate is. 

First, let us assume that one of the two energy transfer mechanisms (r^'''(ijj) or V'^{u!)) is small, 
that is we replace for instance r'^'''(a;) by tT'^'''{uj) with r ^ 1. In this configuration, the energy transfer 
produced by the bottom fluctuations are small compared to the one produced by the surface fluctuations, 
so that 

inf ((-r^-^M + ASM)x,x)j^„<^,, 



T->0 



xes' 



which means that the energy dissipation is caused by the cumulative dissipation coefficients A^'''(ijj) + 
A^'''(w), and only the energy transfer mechanism r'^'^(a;). Moreover, if we assume that the energy transfer 
mechanisms (V'^iuj) and r'='^(a;)) are both small, that is we replace r='"(w) + r'='''(a;) by riV-^iuj)) + 
r'=^''(w)) with T < 1, we have 



lim a: 



.(w) = min A*;'" 
je{i,....N{u)} ^ 



M + AfH. 



According to (76 1, this is the smallest value that Aoo(w) can take. This result is quite natural since the 
coupling process on the energy transfer between propagating modes is small, so that the decay rate of 
the energy of a particular propagating mode j is given by its own radiative coefficient A^''*(tj) + A^'''(a;). 
Then, the decay rate of the total energy carried by the propagating modes is given by the minimum of 
those radiative coefficients. Consequently, if there exists A^^'*(w) + A^^:^^(cj) — 0, we have 

limA^M^O. 

The reason is that the energy of the joth propagating mode stays approximately constant with a weak 
transfer of energy to the others propagating modes, and 



r— !-0 T 



inf (( - {T^'^co) + T^^\u:)))X, > 0, 



where V ^ \ X e S_ 



N{u) 



suppX C { 1 , . . . , iV(w) } \ supp( A^^" (w) + A'='''(a;) ) I , 



since the normalized 



eigenvector (l/^iV(a;), . . . , l/y/N{u})) of r'='''(a;) + r^^''(w) does not belong to V. 

Second, if we assume that one of the two energy transfer mechanisms (r'^''*(a;) or ^'^{iv)) is strong, 
that is we replace for instance r'^''*(w) by ^r'^''*(a;) with t <C 1. In this configuration, the energy transfer 
produced by the surface fluctuations is much stronger than the one produced by the bottom fluctuations, 
so that 

lim A;,(tj) = A(w). 

T— >0 



According to (76 1, it is the largest value that Aoo(w) can take. The strong energy transfer mechanism 
given by just one of the two mechanisms produces a mixing effect, and then a uniform distribution of 
energy over the propagating modes, so that the energy dissipation rate is given by the arithmetic mean 
of the dissipation coefficients. The same result holds when the two energy transfer mechanisms are 
strong, that is we replace r'^'*(a;) + r'^'^(a;) by ;^(r'''^(a;) + T'^'''{uj)) with r ^ 1, but a more convenient 
way to obtain this result is to use a probabilistic representation. In fact, we have 



1=1 
1=1 



exp 



exp 



A 



L 



Lf-r 



A 



N(u,) 



A" 



,M 1 



Yr 



N{ui) 



where (Y^ 



't>0 



is a jump Markov process with state space {I, . . . , N(uj)} and intensity matrix 



r^'''(w) +r'^'''(a;). Using the ergodic properties of (Y^ 



't>0 



we obtain 



lim 'T['\uj, L) 



exp I 



-A(w)L) 



N{uj) 



E 

1=1 



a;,o(w)P, 
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where the first term on the right of the previous equality represents the energy dissipation, while the 
second one represents the energy equipartition. As a result, a strong energy transfer produces the 
strongest energy dissipation that we can expect. 

Finally, in the case of negligible radiation losses, that is we replace A'^'''{uj) + A'^'''(aj) by t{A'^-^{uj) + 
A'='''(cl))) with T <C 1, we have 

limA;,H = 0. 

In fact, if the radiative loss is negligible, the energy transfer mechanism becomes dominant, and we can 
show that 

VL>0, sup ||7^-(w,z)-7;°(c^,z)||2,R«(.) =0(r), 
ze[o,L] 



where T°(w,.) satisfies (73l without the term A^'''(w) + A'^'''(a;). In this situation 



1=1 

so that the total energy is conserved (see Figure [5|, there is no radiation loss, and we have 

lim -A;,H-AH >0. 

T— >0 T 

Moreover, for long propagation distance, as it was already observed in [TTJ [TSl [T71 [1^] , a mixing mecha- 
nism of the propagating modes coupling leads to a uniform distribution of energy over the propagating 
modes, that is 

lim Tf{uj, L) = — — V \Zifl{u 

L-S-+00 ■' N[(jJ) ^-^ 



1=1 



Conclusion 

In this paper we have analyzed wave propagation in a shallow-water acoustic waveguide with a randomly 
perturbed free surface and uneven topography, and the resulting coupling mechanism between the three 



kinds of modes (propagating, radiating, and evanescent). The diffusion system obtained in Section 4.1 
describing the mode coupling mechanism can be used to simulate numerically wave propagation in a 
waveguide with perturbations of its geometry 4 . 

We have shown that the evolution of the forward propagating mode amplitudes can be described in 



term of a diffusion process (Theorem 4.1 Theorem 4.2 and Theorem 4.3 1 taking into account the main 



coupling mechanisms: the coupling with the evanescent modes induces a mode-dependent and frequency- 
dependent phase modulation on the propagating modes, the coupling with the radiating modes, in 
addition to a mode-dependent and frequency-dependent phase modulation, induces a mode-dependent 
and frequency-dependent attenuation on the propagating modes. In other words, the propagating modes 
lose energy through the radiating modes into the bottom of the waveguide and their total energy decays 
exponentially with the size of the perturbed sections. Moreover, we can express the decay rate in terms 



of a variational formula over a finite-dimensional space (Theorem 5.1 1 
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6 Appendix 

6.1 Spectral Decomposition in Unperturbed Waveguides 

This section is devoted to the presentation of the spectral decomposition of the Pekeris operator + 



k {ijj)n {x) , where the index of refraction n{x) is defined by (33 1. The complete spectral analysis of this 
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operator is carried out in |38j . Here, H = L^(0,+oo) is equipped with the inner product defined by 

(■ + 00 

W{hi,h2) e H X H, {hi,h2)^= / hi{x)h2{x)dx. 

Jo 

Definition 6.1 Let us note by R{uj) the Pekeris operator defined by 

R{uj){v)^^y + e{uj)n''{x)y VyeP(i?M), (77) 
which is an unbounded operator on H with domain 

V{R{uj)) = i7o (0, +oo) n iJ^(0, +oo). 
The following result regarding the spectral decomposition of the Pekeris operator is proved in [55] . 

Theorem 6.1 R{uj) is a self-adjoint operator on H . and its spectrum is given by 

Sp{RiLo)) = (-oo,fc2(^)] u{/32^(^)M,...,/32H}, (78) 

where Vj G {l, . . . , N{uj)} , the modal wavenumbers /3j(w) are positive and ordered in an decreasing way 
Moreover, Let be the resolution of the identity associated to R{uj), such that £ H and \/r G M, 



n^(r,+oo)(y)(a;) = ^ {y,(t)j{u:, .)) ^<t)j{uj,x)^r.,+oo) (/3j(w)^) 

(y, (j)^{uj, .))^(l)^{uj, a;)d7l(_oo,fe2(„))(r), 



and Vy G X'(i?(w)), Vr G R, 



n<.(r,+cx))(i?(a;)(y))(x)-^ /3, 

+ / l{y,4'Ji^^,x))JJ4l^{u},x)d'yl(^_ao,k^u))ir). 



Let us describe more precisely the decompositions ( 79 ) and ( 80 1 



Discrete part of the decomposition Vj G {l, . . . , N{ijj)^, the jth eigenvector is given 
(/)j(w,x) 

where 



Aj(uj)sm(aj{LL))x/d) if 0<x<d 

Aj-(w)sin(crj(a;))e-^^(")^ if d<x, 



a,{u;)^d^nlk^u;)~/3^{u;), (^.{u:) ^ d^ Pjiu^f - k^{uj), 



and 



Aj{uj) 



2/d 



\ 1 _i_ si^'^ ("^j ('-^)) s\n{2a j {uj)) 
\ ^ CTM 2<T,(a;) 

According to [38' , cri(w), . . . t<Jn(ui){^) &re the solutions on (Q,nik{uj)d6) of the equation 

tan(y) 



^(nifcd6')2 -y2' 



(79) 



(80) 



by 



(81) 
(82) 



31 



such that < ai{uj) < ■ ■ ■ < ai\j(^^){Lu) < nik{uj)dO, and with 6 = y^l — \/n\ = ^1 — cl/cl- This 
last equation admits exactly one solution over each interval of the form (tt/2 + (j — l)n, 11/2 + j'tt) for 
j G {1, . . . , N{uj)}, so that the number of eigenvectors is 



N{io) = 



nik{uj)d 



(83) 



where [•] stands for the integer part. 



Continuous part of the decomposition For 7 e (—00, fc^((jj)), we have 

{A^{uj) s\n{r]{uj)x/d) 
A,{uj) (smiviu^)) cos (^(w) ^) + fg cos(ry(^)) sin (e(^) ^) 

where 
and 

Let us note that we have 



r]{uj) = d^nlk'^iuj) - 7, ^{uj) = d^k'^{u}) -7, 



if 0<x<d (84) 
if d < X, 

(85) 



tt{S^'^{uj) sin ■^(77(0;)) + rj'^iio) cos2(7y(ijj))) 



AH 



-00 v^|7|i/4- 



(86) 



We can also remark that (j)-y{uj, .) does not belong to H, so that (y, (/)-,(w, .))^ in (79l and (8OI is not 
defined in the classical sense, but in the following way 



M 



y{x)(j)j{uj, x)dx in ( — 00, fc^(w)) . 



As a result, we have 



\\y\\H= 2^ \{y:<Pji^^-))H\ + / \{y:'f>j{^,-))H 



d7, 



and then. 



y 



is an isometry, from H onto = C^^"^^ x L^( — 00, A:^(a;)). 

In this paper the pressure field (12 1 can be decomposed according to the resolution of the identity 
introduced in this section. 



Po{uj,u,v)^ ^ p.j{uj,u)(j)j{uj,v)+ / p^{uj,u)(j)^{uj,v)d'y+ / Pj{u;,u)(j)^{uj,v)d-f, 

Jo J-00 



radiating modes 



propagating modes 

where the amplitudes of the propagating modes are defined by 



evanescent modes 



Vj e {1, . . • ,iV(w)}, Pj(cJ,u) = n„({/3j(w)})(po(w,u, •)) = / po{uj,u,v)(j}y{uj,v)dv, (87) 



+00 





and the amplitudes of the radiating and evanescent modes are defined by 

V7 € (— (X), fc^((jj)), pj{uj,u) = lim / po{uj,u,v)(f>j{uj,v)dv in ( — cx), fc^(a;)) 

Consequently, the mode amplitudes are the projections of the pressure field po{uj, u, v) over the eigenele 
ments of the Pekeris operator defined by ( 77 1 . 
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6.2 Proof of Proposition |2.2 



As we have remarked right after the statement of Proposition |2.2| the conformal transformation (f> can 
be extended as a homeomorphism from the closure of Dq onto the closure of Dy, ■ Consequently, letting 
w — >■ in the definition of z{u,v) in (16 1, we obtain the following equation 

which is well defined since V/ has a compact support and z{u,0) goes to +oo (reps. — oo) as u goes to 
+0O (reps. — oo), and this equation also can be recast as follows 



q{u) - Q{q)(u) = q{u) 
where q{u) = z(u,0) — u, q{u) = —U{Vg){u), and 

1 



(89) 



Q{q){u) 



-p.v. 



F/(u + g(7i))-F/(u) ^,^ 



The goal is to invert (89 1 in _L?'(IR) (p > 2) by showing that Id — Q is homeomorphism, where Id stands 
for the identity map. However to show that we need to know first ii q G LP(IR). In fact, Vu > one can 
show that 



|($ - Id){u + iv)\Pdu < +00, 
where <1? is the conformal transform, so that according to [37j Theorem 93] we have q E LP{M.). Moreover, 
\\QmLpm<V~^K\\q\\L.m and Wqh^m < e'^^-'^^K 



thanks to the boundedness of the Hilbert transform operator in 



so that ||g||LP(R) < K thanks to 



(89 1 for e small enough. Here, K denotes a deterministic constant which is independent of e. Now, let 
us consider Q the application from LP{M.) to itself defined by 

' 

— n times 

This application is well defined since LP{M.) is a Banach space and ||Q(OIIl!'(k) 5; v^^''^II^IIl!'(r)i and we 
have 

(/d - Q) o e = Id and 6 o (/d - 6) = Id. 

Moreover, one can easily check that both {Id — Q) and O are Lipschitz function, so that {Id — Q) is an 
homeomorphism. As a result, we obtain from (|89| that 



q{u) = e{q){u) = q{u) + Q{e{q)){u) with \\Q{e{q))\\LH 
which concludes the proof of this proposition. 



< e 



i-i/p 



6.3 Proof of Proposition 2.3 



First, under Assumption ^ and thanks to the Holder's inequality, we have 
\x{u, v) — v\ < ( — 

\ TT 



■V f \x{u,0) — x{u,0)\ , , 

' ^idu+|a;(M,0)| 1(„<^) 



{u — u)^ + 
x{u,0)\Pdu 



{{U - m)2 + w2)g 



du 



1/9 



'-(i'>m) 



<,i/2-i/PcJl 
Vtt 



du 



\u\>l 



\u\<l/v ^ 



1/2-1/p C2 
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since p>2, sup„>o \U{V^'){u)\ + \U{{d/du)V^){u)\ < Ce^/^-i/p, and 



\z{u,v) - til < (- 



|u — u||a;(u, 0)1 1 f \u\\x{u + u,0) — x{u — u,0)\ 



\u-u\>i {u-u)^ + v^ 



dii -\ 



\x{u,0)\Pdu 



i/p 



|«|>i 



u — u 

1/9 



|u|<l 

{"du 



du]l 



1/9 



L(D>/i) 



\u\<l/v ^ 



For the second part of the proof, V/i > 0, we have 



+ CSO 



+ 00 



E[|z(u,w) - < 2e / dv I du I duR^{u - u) 



fsieu)fs{eu){u ~u){u- u) 



L/e 



+00 



< 2e / dwRs{w) I dv- 



^ l2/9 



It 



+ t;2 



+00 



{t2+ 1 



du 



2/p 



+ e2(i-i/9)c 

f^~|~00 /'~]~CX2 1 T~ 

dwRsiw) dv-ln(^-^ + l^ +e^^^-^/'i'>C, 



+00 



dv- 



1 



1-2/p 



< e 



with p > 2 and 1 < q < 2, and i?^ is the autocorrelation function of the surface fluctuations Vg introduced 
in ([t]). However, 



+00 



dv — In 



L 



+ 00 



1 , r L 



V \v 



1 , / i 



1 <e/ dv-ln{ — + l]+e dv-ln { — + 1) < C{e + e'-") 



en 



V \V 



for some i' G (0, 1). Now, 

^+00 



+ C30 



E[\x{u,v) ~ v\^\ < 2e / dv du duRs{u~u) 



fs{eu)fs{eu)v'^ 



-\--oc 



dvv^E 



2n 



(u — {t)2 + ?j2 
+00 rL/e 



/■L/e r+oc l-L/e ^ 

< 2e / dwRs{w) I dv I du^ 

Jo Jr, Jw U^+V 



+ e2(i-i/9)c 

r'+oo 



L i2/g 



+00 



-du 



2/p 



+ 1)P 

<2e/ dwRs(w) dw-arctanf— +e2(i-i/«)C, 



+00 2^ 



yl-2/p' 



with p > 2 and 1 < q < 2. However, 

^+00 



dw-arctan^ — ^ —^J du- arctan ^ + e J dw-arctan^— ^ < C^e + eln ^-^^ , 
which concludes the proof of Proposition |2.3| thanks to the first part of the proof. 



6.4 Proof of Proposition 3.2 



The proof given in this section follows closely the one proposed in [17, Section 4.3], but we still present 
the main ideas. First of all, we recall that O^j o Hi^(— 1/^, — ^) (^pq(uj, ., m)) represents the evanescent part 
of the pressure field po{u>, ., u), where Q^j and are defined in Section [6. 1[ In this section we consider 
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the Banach space F = 1/^, — ^) equipped with the norm = J_i/^\y-y\d"f- According to (39l 

and the radiation condition given by Assumption [3] we have 



(90) 



where g^(u>^u) is defined by (40), and 



|^|l/4 ' 



|7|«1 



■(-1/? 



-?)(7)- 



Therefore, substituting (40) into (90), we obtain the following relation in (C([0, +oo), F) , ||.||oo,f) 

(M-<i>")(e„on^(-i/e,-o(po(^, .,.))) =p(c^,.) + pi(^,-), (91) 



where ||y||oo,F — sup2>Q ||y(z)||i?, Vy G C([0, +oo), i^). In (91), is a hnear bounded operator, from 
(C([0,+cx)),F), ||.||oo,f) to itself, defined by 



fc2(w) 



2Vl7l ^ ^-i/C 



,(a;, u + i;)2/y (w + w)d7'e"V l^ll^'ld^; 



and 



1=1 



7' 



,1/4 



(ay(a;,u + z;)eVT"("+«) + Sy (a;, w + w)e-VT'("+''))d7' e-vl^H^ldu 



Vm e [0, +oo) and for almost every 7 e (— ^0- We can check that the norm of the operator is 
bounded by 

g3/2-l/p J 

- ( sup I K (u) I + sup I — (w) I + sup I H (u) I) , 



so that lime_j.o V{Id — (f>" is invertible) = 1, and the condition {Id — <i>" is invertible) is satisfied in the 
asymptotic e — ?> 0. As a result, on the event {Id — ^/e<f>'^ is invertible), we have 



e„on„(-i/e,-0(po(w,.,.)) = (w-*") \p{lo,.) + p,{co,.)) 

= p{uj,.)+pi{uj,.) + <i>^{pi{uj,.)) 

+00 

+ J2{^-yipi^,-) + '^^iPii^,-))), 



so that 



o n„(-i/e, -0{po{^, ■)) - p(t^, ■)+Pi{^, •) + 'i>"(Pi(c^, •)) 



O {sup\Vs{u)\ + sup\—Vs{u)\+sup\Vb{u)\) sup ||a(a;, M)||«je + \\b{uj,u)\\n^ 



^ u u dtL 

in C([0, +00), F) . Moreover, introducing 



(92) 



;^.,2(a;,.)^fM ^^'(";^+-) (^,(^,,)e.M...)+SKo.,u)e-'A("+")) 
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Vm e [0,+oo), wc have according to (36 1, (37), and (38), 

e3/2-i/p d s2 

\\p{uJ,.)-p2{l^,.)\\oo^F < C (sup|V;(-u)| +SUp|— V;(u)| +SUp|yb(M)|) 

X ( sup MLJ,z)\\n- + \\H^,z)\\h- + \\e^ oU^{-l/(,^^){po{ij, .))\\f). 



Therefore, since Vg, Vb and their derivative are assumed to be bounded, and according to (92), we obtain 

o n^(-i/e, -0(po(^, .)) - P2(c^, .) + Piiuj, .) + <i>"(pi(w, .)) 

+ 0[ (sup|y,(w)| +sup| — +sup|yb(u)|) sup \\a{oj,z)\\n^ + \\b{aj,z)\\-H^ 



3.2 



in C([0, +00), F) , which conclude the proof of Proposition 
6.5 Proof of Theorem 14.11 

The proof of this theorem is based on a martingale approach using the perturbed-test-function method 
and follows the method used in [IT]. However, the process (a^''^(M))^^p is not adapted with respect to 
the filtration JF^^ = where Tu is defined by ([s]). In fact, neither the random operator nor 

Gj"^ are adapted with respect to the filtration J-^. The proof of this theorem is in three parts. The 
first part of the proof consists in approximating the random coefficients H^'| and G"^ with new ones 
from which we can use the mixing property. The second part follows the ideas of .22j and consists in 
simplifying the coupled mode equation and introducing a new process for which the martingale approach 
can be used. Finally, the third part of proof of this theorem is based on a martingale approach using 
the perturbed-test-function method and follows the ideas developed in [5l [T7] . 

6.5.1 Approximation of the Random Coefficients 



This section is devoted to the approximation of the random operators H^, ^ and G^ ^ defined by (45 1, 
(57 1, (58), and G^ ^ ^y = for 7 G (5, fc^(a;)). In fact, these operators do not satisfy mixing properties 
coming from the random perturbations of the waveguide surface. The conformal map (16 1 involves the 
complete trajectory of the random surface of the waveguide. 
Let ao e (1/2,1), 



Vsiu/e + w)fs{u + ew) 

'\w\<l/e''0 

2 f Vs{u/e + w)fs{u + €w) 



1 /9 , /u\ „ , , N , / ,N , / ,N " / Vs(u e + w)js(u + ew) , 

1/2,. _ ^2^^ _ d _ / ujy 



f +00 



Mv)Uv) / o : „ -dy^dv (93) 



and 



Ps.t = (t{C^'^''{u), s<u<t), (94) 

the cr-algebra generated by C^/^'*^. Let us remark that Pg.t is still ^^-mixing. In fact, for s > 2/e"°, if 
A G Tt+s,+oo and B e -Fo.t, then A G .Ft+s_ei-c«o.+oo and B e ~^o,t+ei-"o, and therefore 

sup \V{A\B) - P{A)\ < <j>{s - 2/e"°) = ^e(s). 

t>0 

A&ft + s, + o^ 

According to we have the following results regarding mixing processes : 

|E[Ci/2-^(i + si)|j£-o,t]| <2^,(si), (95) 

and 

\E[C'/^-'{t + ,si + s2)C'/^'%t + s,)\^o.t] - E[Ci/2.^(t + + s2)Ci/2'^(t + s^)] \ 



(96) 
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Vsi,S2 > 2/e"». 

Let us also consider 

2{ni - l)(l)r{d)(t)s{d) 



dxP f Vs{u)fs{eu) 



{u-uY + (P"^J {{u~uY + (F) 



-du 



-du 



3d / 



{c^r{d)^'s{d) + cl>'M)Ud)) 



-n ^du 

n J (u — m)2 + 



d r U{Vs{-)fs{e-)){u)j^Vs{u)fs{eu) ^^ 



+ 00 



{u~u)Vs{u)fs{u) 

{{u — u)2 + 1)2)2 



z^(y,(.)/.(e-))(ii)i^v;(zi)/,(«i)^^ 



4i;2 fUiV.i-)fsie-)){u)^^V,iu)f,ieu)^^ 2v^ f Vs{u)fs{eu) 



((u — u)2 + u2)2 



{u ~ '(i)2 + ti2 



-du 



{U — U)2 + 1)2 

V;({l)/,(ei2) 

((ft — u)2 + ti2)2 



du 



Vs{u)fs{eu) ^ x2 ^ 4^4 /■ F,(U)/,(6U) ^, 

(u-m)2 + ?;2 7r2 V 7 ((u - u)2 + w2-)2 " 



(97) 



In the expression of C^/^'*^ and C^'*^, that is (93) and (97), all the terms involving a (f>rid) produce, 
through the conformal transformation (16), a mode coupling at the bottom of the waveguide, while all 



the others terms produce a mode coupling all along the transverse section of the waveguide . As a result, 
the random perturbations of the waveguide surface produce, through the conformal transformation, 
perturbation in the body and at the bottom of the waveguide. 



Lemma 6.1 V(r, s) G {1, . . . , A^} U (C, fc^), 



have 



C;,(u) = CU^'^iu) + e Clfiu) + o(e). 



uniformly with respect to u ^ (0,L/e), where C'^{u) is defined by (32) 



This Lemma gives an approximation of the random operators H°'| and G"^ defined by (45), (57), (58 1, 
and G"'^ ^ = for 7 e (^, fc^(w)), by new operators H^'^ and G^'^ defined by 



H„.,«(")(y) 



yie 



2 N ^1/2,6/ X 



N 

E 



71 TTtA 

cj/^^(u)c;/f''(u+u) 



Vie 



i/3i (m+m) - •y' 1 7' 1 1 " I - i/3j " 



ik^ 



N ri\.t 



[E 



where C^l'^^^u) = if u ^ [0,i/e] and g!^",^ = for 7 e (C,fc2(cj)). The following Lemma is a direct 
consequence of Lemma |6.1| 

Lemma 6.2 We have 

uniformly in u £ [0, L/e] and for each y e "H^. 
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This second Lemma will be used in the next section to approximate the forward mode amplitudes a'^'^ 
unique solution of Equation ( 56 I . 



Proof (of Lemma 6.1 ) The proof consists only in doing expansions of the perturbations produced by 



the conformal transformation. 

^+00 



C^u) 







where Jg is defined by ( [28| ), and 



+ 00 



(l)r{v)(l)s{v){J^{u,v) - l)n'^{v)dv, 



-\-oo 



(t>r{v)(t>s{v){nl{u,v) -n^{v))dv, 



(j)r{v)cj)s{v){Je{u,v) - l){nl{u,v) - n^{v))dv. 



First, 



Second, 



where 



+00 



(j)r{v)(j)s{v)rL^{v) 



(i)r{v)(j)s{v)n^{v) -^(^ 



4v^ I Vs{u)fs{u) 2 r Vs{u)fsiu) 

f f {u-u)Vsiu)fsiu)^A^ 



du 



2 fU{Vs{-)fs{e-mu)£Vs{u)fs{eu)^^ , Av^ fU{Vs{-)fs{e-)){u)j-V,{u)Ueu)^^ 



dz 



2t;2 r Vs{u)fs{eu) 



du 



-du 



Vs{u)fs{eu) 



((u — It)^ + u^)^ 



-du 



J (u — uy + J ((u — u^ -\- 

1 



du 



TT 

0(e). 



{u — w)^ + 



{{u — w)2 + ti2)2 



{1 - ni)(l)r{v)(t)s{v)dv -\- / {ni - l)(t)r{v)(t)s{v)dv, 



^i,e = {v < d, x{u,v) > c?} and A2.e — {v > d, x{u,v) < d}. 
The real part x{u, •) of the conformal map (|16p is a nondecreasing bijection since 



sup \dvx{u,v) — 1| < ^/eK. 

u>Q 

ve[n-M] 



Moreover, 



\x ^{u,d) — d — {d — x{u,d))dyX ^{u,x{u,d))\ < -\d ~ x{u,d)\'^ sup {d^x ^{u,v)\, 



with 



Then, 



sup sup \d'^x ^{u,v)\ < K^/e, and sup\d — x{u,d)\ < K^/e. 



u>0 ve[d,x{u,d)] 



u>0 



{1 - ni)4>r{v)(l)s{v)dv + / {ni - l)(j>r{v)(l>s{v)dv + o{e), 
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where 
and 

Y'{u,d) 
Consequently, 



A2,e^{v>d, v<d + Y%u,d)}, 



x{u, 0) 



rdu 



1 

1 - - 



TT J ({t — + rf2 L nj (m — + rf2 TT J (({t — + cJ^) 



0) 



Z?2.^(u) = ^/e 2(711 - d)- 



du 



+ e2(ni - l)(/)r('^, d)4>s{uJ, d) 



d r Vs{u)Meu) 

TT J (m - m)2 + d2 

r d3 /• Vsiu)fsieu) 



TT^ J {u- uY + d'^"'"' J {{u-uY + d^) 



du 



rdu 



d ( f VMfsieu) , x2 d fU{Vs{-)fs{e-)){u)£V.{u)Ueu)^^ 



7r2 V 



(u — u)2 + d^ 



du 



+ (M^, d)^',{uj, d) + d)) 



[u — u)2 + ^2 

t4(ii)/«(eu) ^ x2 



du^ +o(e). 



Finally, according to the previous step 

r4d2 



TT y ("S — u)^ + d"^ 

Vsiu)fsieu) ^id f V,iu)Ueu) 



du 



TT J ((u-u)2 + d2)2 T:J {u-uY+d^ J TT J (u-u)2 + d2 



c?u + o(e). 



To finish the proof of this lemma, we have to approximate in D^''^{u) and D^''^{u) the terms of order ^/e 
to obtain a more convenient form to exploit the mixing properties of the random perturbations. 
Let ao G (1/2,1). First, 



Vs{u)fs{eu) 



du — 



and then 
In the same way, 

f + 00 



|-!«|>l/e°o 

M + tU>0 



{u — m)2 + d"^ 
Vs{u + w)fs{e{u + w)) 

w2 + ^2 



Vs{u + w)Ss{e{u + w)) 



dw 



< K 



+ 00 



dw. 



-dw < Ke"°. 



Vs{u)fs{eu) 
(u — lt)2 + i;2 



dudv 



+ 00 



4ir{v)(j)s{v) 



l/(de°o) W2 + 1 

Vs{u + w)fs{e{u + w)) 



yj2 ^2 



dwdv. 



so that 



+ C30 



lu)|>l/e°o 
u+u;>0 



< is: 



< 



Vs{u + w)fs{e{u + w)) 
|0,.(^;)0.(i;)| '■+°° 





+ 00 



l/(De°o) 



dwdv 
1 

0) ^2 + 1 



dwdv 



\Mv)My)\dvKe° 



Let us note that f^°°\(f>r{v)(l)s{v)\dv is well defined thanks to Assumption [s] and to the form of given 
in Section |6.1[ The last term can be treated exactly as the previous one. That concludes the proof of 
Lemma [6A1 □ 
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6.5.2 Approximation of the Transfer Operator 



The forward mode amplitudes a'^'^ unique solution of (56 1 cannot be studied directly using the perturbed 
test function method introduced below because the random operator H^'^ and Gg"^ are not adapted to 
the filtration of the random perturbations. Consequently, we approximate a*^'^ by a new process from 
which we can exploit the mixing properties of the random perturbations. 
Let us introduce d'^'^ the unique solution of the differential equation 

^d^Hu) - ^Hr, (^) d^Hu) + (G"")(u)a^-«(«), (98) 
with a<='«(0) = 4 defined by ([24]), and where {G"") is defined, \/y G H^, by 



/"'!^ r G^.;)(.)cos(/3,.)e-^^^...,,, + ^Gf 



Vj e {1, . . . , iV}, and (C"") (y) = for almost every 7 e {£,, fc^). Here, 



r<Wf \ A( 1^2J,2^^^ ,2/ nO'^ [[ Rs{z + [w - w')d) , 



^ // <Pj{v)'f)-i{v)ci)j{v')(i)^{v') [ [ /^f^ ^ . ^ ^ j dwdw'dvdv' 

J J J J {W^ + V'')(W' +V' ) 

^ ^ (i)j{v)cj)^{v)cj)j{v')(i)^{v')v'^v''^ 11 ^''^^ ^ „ ^ ^ I dwdw'dvdv' 



—(t>jid)(t>^{d) I (j)^{v)(j)^{v) I ,^2^f ,2 ^ dwdw'dv 

TT J J (w^ + w^)(w' + a ) 

16(ni — 1) , / ,N , / ,N f I f \i ( \ 2 f Rs{z + WW — w'd) 

TT 



-(j)j{d)(j)^{d) / <?!)j(w)0^(w)w / , , , 2^2/ ,2 , dwdw'dv 



<^.(-)'^.(-)'A.(-)<^.(-)« yy (^2 + ,2)2(^,2 ^,,2) 



and 



^ // ^ dwdw' ^ // ^ dwdw' +U(R,)(0) 

TT^ JJ iw^ + d^){w'^ + d^y n^JJ (ii;2 + d2)(u,'2 + d2) ^ \ 



, / ^ , / N 2/ ^^4^'^ ff ww'RJw-w') , , , 2,,,-- 
+ / 0, w 0, w w // ,2 , 2W 2 , ' dwdw' + -UiR.m 

2t;2 ff R,{w-w') 

+ -(;2)(^i;/2 _|_ y2y 

Rs(w-w') , , , 4^4 ff RJw-w') , , , 

^ ' -dwdw' + ^ // ^ dwdw' 



TT^ 77 (^2 + w2)(-y;,2 _|_ ^2) 7r2 77 (-y;2 _|_ ^2 )2 (-^,,2 ^2 -)2 

where i?^ is the autocorrelation function of the random process , and Rg is the autocorrelation function 
of the processes Vg and ^ Vg , that is j^Rs- 

We have the following proposition that describes the relation between the two transfer processes 
o'^'^(m) and d'^'^{u). 



Proposition 6.1 



V77>0, limP sup \\a''^{u) -d''^{u)\\l^^ > r]] ^0. 

\uelo,L] « / 
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Thanks to Proposition 6.1 one can study the new process (^a'^'^u)'^ instead of (a'^'^u)) Let 



us remark that the new process (a'^'^(u)) , is adapted to the fihration J-q ^ defined by (94 1 and 



\\&'Hu)\\ 



l^ji Vu > 0. Consequently, {a'^'^ {u)) can be studied using the perturbed test 



function method. 



Proof (of Proposition 6.1) This proof is in two steps. The first step consists in using Lemma 6.2 



to 



approximate the random operator ^ and ^ by ^ and ^ from which we can exhibit the mixing 
properties of the random perturbations. The second step consists in using the ideas developed in 
to exploit the fast phase of the random operator G"*^ . The proof of the second step follows closely the 
proof of Proposition 6.3 in [T5]. ■ 



6.5.3 The Perturbed Test Function Method 

In this section we study the limit in distribution of the process 



[of^'^ u>o unique solution of the 
differential equation (98), using the perturbed-test-function method and a martingale technique. The 
proof is in two steps. First, we prove the tightness of the process {a'^'^ {uf^ afterward we identify 
all the subsequence limits thanks to a well posed martingale problem. According to Theorem 4], 
the proof of the tightness is decomposed in Lemma |6.3| Lemma |6.4|, an d Lemma |6.5[ The subsequence 
limits are identify, first by a complex martingale problem (Lemma |6.6[), and then by a real martingale 



problem (Lemma 6.7) to use the classical uniqueness results in |3D] for instance. 



Now, let us defined the good space to apply the asymptotic analysis. Let r 



the closed ball with radius r, and n > 1} a dense subset of Br^-Hf ■ We equip Br^Ut; with the distance 
dBr H defined by 



+ 00 ^ 

|(A-/^,5n) 



V(A, /i) S Br.n^'^- As a result, {B-^^ , dg^ „^ ) is a compact metric space. From the definition of the metric 
d-Br H ' w2 have the following criterion. 



if and only if 



Theorem 6.2 A family of processes (X'^)£g(o,i) tight in C([0, +oo), (S^. 
{{X%X)^^)^^f^Qj^j is tight onC([0,+oo),C) VA e H^. 

This last theorem looks like the tightness criterion of Mitoma an d Fouque [301 HH]- For any A G 

the family (a'^'^(.))c is tight on 



(^d'^'^{u), X)^ . According to Theorem 



6.2 



■Hj, we set a^'^(ii) 

C([0,+cx)),(6^,«j,de,„J) if and only if the family («!'*(. ))e is tight on C([0,+oo),C) VA £ U^. Fur- 
thermore, (a'^'^(.))e is a family of continuous processes and then it is sufficient to prove that, VA G H^, 
(a^'^(.))e is tight in the space of cad-lag functions I'([0, +oo), C) equipped with the Skorokhod topology. 
According to proposition [6A] to study the convergence in distribution of the process (a'^'^(.))^ it suffices 
to study this convergence for (a'^'^(.))^. 

We consider the complex case for more convenient manipulations. Using the notion of a pseudogen- 
erator, we prove tightness and characterize all subsequence limits. 

We recall the techniques developed by Kurtz and Kushner. Let A^*^ be the set of all J^'^-measurable 
functions f{u) for which sup„<^ E [|/(u)|] < +oo and where L > is fixed. The p — lim and the 
pseudogenerator are defined as follows. Let / and /"^ in A^'^ V(5 > 0. We say that f = p — Vang f^ if 

supE[|/'^(u)|] < -I-C50 and lim E[|/*(u) - /(w)|] = Vu. 



The domain of A" is denoted by 2? (^'). We say that f eV {A") and A' f 
and 

p. -K[f{u + S)]-f{u) 



g if f and g are in T> {A'^) 



lim 

5^0 



.9(") 



0, 



where E^ is the conditional expectation given J-^ and J-^ 
the following theorem. 



J-u/e- A useful result about A"^ is given by 
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Theorem 6.3 Let f e ViA"). Then, 

M}{u) = f{u)~ r A'f{v)dv 
Jo 

is an {J-'^)- martin gale. 

In what follows, we consider the classical complex derivative with the following notation: liv — a+ij3, 
then dy = \ {da — idp) and — [da + idp). 

Proposition 6.2 VA G Ti.^, the family (a^'^(-))£g(o i) tight in I? ([0, +00), C). 
Proof First, we easily obtain VT > 0, 



lim lim ] 



M-i.+oo e-s-0 \o<„<T 



sup |al'«(u)| >M) = 0, 



since {a'^'^{))e is a bounded process. To show the tightness of the process (a'^'^(it))^^p, according to 
[2S1 Theorem 4], we need to show the three following lemmas. Let A G T-L^, / be a smooth function, and 
/o^(w) = /(a^/(u)). We have. 



where 



H 



and G^{¥^^{u)) = {{G^^){u)a^^^{u),\)^^^ 
mu) = -^9./ (al-«(«)) 



Let 



+ 00 



+00 



iI.(a^'«(u),GV2'^(^), 
i/A(a^^«H,Gi/2,.(^) 

Here, E^ stands for the conditional expectation with respect to the filtration Jq defined by (|94| 



Lemma 6.3 VT > 0, linig supq<„<7.|/J^(u)| = almost surely, and sup,„>Q E = 0{e^^^ ^) 

V77 > 0. 



Proof (of Lemma 6.3) Using a change of variable and thanks to (95), we obtain 



\n{u)\<K^e , 
Jo 



(j)(w)dw 



-\-OQ 



V-HCl + l/{2w))dw 



-j-oo 



+ A'yiln(l + l/id^e^"")) + K^e / ^{w){\n{\ + \/{w^e^°">)) + l/(2u;))(l + \/w)dw, 



where, thanks to Assumption [s] il> is a function coming from the propagating modes (t>j such that the 
integral in w is well defined. More precisely, ijj{w) ^ uP' as w goes to and decay exponentially at 
infinity. The second and third part of the previous expression come from the fact that in ( 95 1 s has to 
be greater than 2/e"°. Indexes s smaller than 2/e"" are not covered by the mixing property and then 
lead us to explicit computation using a primitive of arctan(?«) given by i/; arctan(i/7) — l/21n(l + w^). 
That concludes the proof of Lemma[6.3|D 
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A computation gives us 
^^(/o + /0(«) = 



where 



, — , — \- w ] aw 
e e / 



H 



+ GA(a^'«(u))" 



/■+00 

^0 



dw 



w I dw 



0(1), 



+ 00 - 



N N 



V / 



(ft-ft)5+i(V7"-ft)srj,_^„^y/ 



4 



^7"'' p»(ft-v/7)S+i(ft'-ft)sT;, 



E / 



iff (T,C,s,S) 



AT Af 



-/3,-)5+i( V^-/3^., )srj,^rj,_^, 



+ / ^ 



+ 



C. 
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ijf (T,C,s,5) 



^/ / — ^' 



„i(/3i-ft)s-i(/3,/-^7i)srj, rjn 



« ''=1 V/3^-\/7W2/3i' 

E [ E 



g«(A-y7r)«-«(/3i'-/9j')s'P;T;7 



E 



1=1 

4 



a 



71^72 



g*(A-y7r)^--«(v^-/3,0^T,T^^d7^ 



JiPl-VrL)S-i{Pl'-Vr2)sr^ rj^ 



Let us note that a quick computation shows that A'^if^ + fi){u) is not necessarily uniformly inte- 
grable. In fact, we have just obtained for E[|^'(/q + the same bound as the one obtained for 

fi{u) without ^/e, so that this bound blows up as e goes to 0. The problem comes from the indexes not 
covered by the mixing property (95 1. To show the tightness of the process d"^'^ we have to introduce a 
second perturbed test function to correct the problem. This second perturbed test function /| is given 
by 
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where 



+00 



Fx{T,C,w)=dJ{Tx) I h[^\t,C{w) (E) C{w + w),w,w + w)dw 



+ 00 - 



+ dvf{Tx) / H^^\T,C{w)®C{w + w),w,w + w)dw 



+00 



+ d'tf{Tx) I H^^\t,C{w) ® C{w + w),w,w + w)dw 
Jo 

r+oo 

+ dlf{Tx) / H'^^\T,C{w)(^C{w + w),w,w + w)dw 
Jo 

+ dwdvfiTx) / H^^\t,C{w) (g>C{w + w),w,w + w)dw 
Jo 



+ dvdvfiTx) / H'^'''{T,C{w)(g>C{w + w),w,w + w)dw. 
Jo 

Lemma 6.4 VT > 0. liirij supg<„<2n|/|(u)| = almost surely, anrf sup„>Q E[|/|(u) |] = 0{e^~'^), Vry > 0. 



Proof (of Lemma 6.4) Using a change of variable and thanks to (96), we obtain 

2 



\f2{u)\<Ke 



+00 



+ 00 



^^/^Hrfw / V(w)(l + l/(2w))dw 



+ ii:e(^ln(l + 1/(^2 e2"o))j + K\ ilj{w){\n{l + 1/ {w'^e^°'°)) + l/{2w)){l + l/w)dwj , 



where ijj has been described in the proof of Lemma 6.3 We recall that indexes s smaller than 2/e"° are 
not covered by the mixing property. These indexes lead us to explicit computations using a primitive 
of arctan(u') given by warctan(i(j) — l/21n(l + w^). That concludes the proof of Lemma 6.4 I 



□ 



Finally, the tightness of d'^'^ is given by the following lemma. 



Lemma 6.5 {A''{fQ + fl + /ijC"),^ G (0, 1),0 < u < T} is uniformly integrable VT > 0. More parti 
ularly, 

^1/o + /r + /2)M = 

9„/(a^'«H)[ / ^ H^^Hd''^{u),E[C^/^^'{0)(g)C^/^''{w)],-,-+w)dw 
+ Gx{d''Hu)) 

+ 00 



G. 



+ 0(1), 



+ 00 



H 



(a^'«(u),E[Ci/2^^(0)®C7i/''^(u;)], 



+ 00 - 



H 



(2) ~ 



dw 
dw 



+00 



H 



- 



H 



(a^.€(w),E[Ci/2,e(o) ® C7i/2,e(u;)]^ 7' 7 + 
(a^'«(u),E[Ci/2-'(0) ® Ci/2^^(zi;)], ^, ^ + , 
(a''«(u),E[Ci/2,^(0) (g, Ci/2.^(u;)], ^, " + i«) 



dw 



Proof (of Lemma |6.5[ ) This lemma is a consequence of a long but straightforward computation. The 
terms of order one in A'^iy f^ + /i + /|) are uniformly integrable. Moreover, the terms bring by /| in 
0(1) can be bounded using the bound obtain in the proof of Lemma 6.4 with ^/e instead of e, so that 
the unbounded logarithm which was a problem in A'^{f^ + fl) can be killed. □ 
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That completes the proof of Proposition |6.2| 



Now, we characterize all subsequence limits (Proposition 6.3 1, using first by a complex martingale 
problem (Lemma 6.6I, and then a real martingale problem (Lemma 6.7) to use the classical uniqueness 
results in [ID] for instance. 

Proposition 6.3 All the subsequence limits of {a'^'^{-))ee(o. i) '^^s solution of a well posed martingale 
problem associated to the infinitesimal generator defined by (59 1. 



Proof (of Proposition 6.3) To do that, we consider a converging subsequence of (a'^'^(.))eg(o,i) which 
converges to a hmit a^(.). For the sake of simphcity we denote by (a^'^(.))cg(o,i) such a subsequence. 

To exhibit the well-posed martingale problem we need to extract all the fast oscillating phase in 
A*^ (/q + /i + /l) which will not be canceled in the asymptotic e goes to 0. 

Lemma 6.6 VA G "H^ and V/ smooth test function, 

f{di{u)) - / a„/(ai(w)) {j^{Z^{w)),\) +(hf{ai{w)){J^{Zi{w)),\)^^ 



+dlf{Zi{w)) (/f(a«(H)(A),A) +a|/(Si(u;))(if(a«(H)(A),A) 



+d^dj{Zi{w)){L{a^{w)){X),\)^ + d,chjf{Zi{w)){L{ai{w)){X),\) dw 



is a martingale, where 



pi. 

J3 



A 



2 



-\T. r^VT.TAz - \ E (r;r + j) T,T,Az, 

1=1 1=1 

1 ^ _ 1 ^ _ 

L(T)(A), ^-Y^T'^TiliXi + -5]p;fT,T,A„ 



(99) 



1=1 



and J^(T).y = K{T){X)-y = i(T)(A)^ = for almost every 7 G i^,k'^), and for (T, A) G Jij. Here, P 



l,s 



pc,s^ ^c,s,i^ ^s,s,i^ ^s,i defined in Section 4.1. 1 



To prove this result, we assume that the following nondegeneracy condition holds. The wave numbers 
j3j are distinct along with their sums and differences. This assumption is also considered in pTHT^f^O] . 
The proof of Proposition 6.6 follows exactly the proof of Proposition 2.9 in pT^ . 

In order to prove uniqueness, we decompose a^(.) into real and imaginary parts. Then, let us consider 
the new process 



Y«(w) 



.Y''^(m). 



This new process takes its values in C/^ x C/c, where t/^ = 
operator 



, where Y^'«(u) = Re{a^{u)) and Y^'^(u) = Im{a^{u)). 

^ X L'^{{^,k'^),R), and we introduce the 





Si — 










r|-i2 




r|-i2 







Lemma 6.7 V/ G 0^(0^ x g^). 



Mj{u) = f{YHu)) - / L^f{YHw))d^ 



is a continuous martingale, where V(Y, A) G {G^ x Q^)'' 



L«/(Y) = -trace {A{Y)D'f{Y)) + (B^ {Y) , D f (Y)) 



(100) 
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with A{Y){X) = Ai(Y)(A) + A2{Y){X) + A3(Y)(A). Moreover, for j e {I, . . . , N} , 



' 3 J 

2 



A 



N 



A 



T,(Y) 



A,{Y){\), = T,{Y)Y^T]f[T]{Y)\] + T2(Y)Af] 

1=1 

N N 

^2(Y)(A), = -Y.Y^V^^fiYlXl+Y^Xj] +T,(Y)5]r;f [Ti(Y)Ai + Tf(Y)A2] 



1=1 



1=1 



N 



^3(Y)(A), = A,^r;-[(Yi)2 + (Yf)2], 



1=1 

and B«(Y) = A^(Y)(A) = A^(Y)(A) = A^(Y)(A) = for almost every 7 G (C,fc2). Moreo?;er, i/ie 
martingale problem associated to the generator is well-posed. 

Proof (of Lemma 6.7) Following the proof of Theorem 4.1.4 in [JD], to prove that Alj is a martingale 
it suffices to show that 

J 

is a continuous martingale with quadratic variation 

< Ml > (u) = £ {A(YHw))iX), X)g^^g^dw. 

Moreover, V(Y, A) e {G^ x C/^)^, we have {A(Y){X), X)g^^g^ > and trace{A(Y)) < +00. According to 
Theorem 3.2.2 and 4.4.1 in [40 , the martingale problem associated to is therefore well-posed. □ 



Finally, to recover rigorously the form of the generator (64 1 for the complex diffusion process a^(a;, •), 



we use formula (100) where 
with 



g(Y) = /(T) and T = Y^ + iY^ 



9Yig(Y) = (9t, + V)/(T) and d^^giY) = i{dT, - d^)f{T). 
In fact, for the drift term, 



(i?«(Y),Ag(Y)) 



N 

E 



N 



r 



c,s 

n_ 
2 



' n 
2 



a;-' 



Y?aY^g(Y) 



r 



2 



Y,^9y^5(Y) 

3 

Y,^9y^5(Y) 



A' 



9 



jOt- + C.C, 



where c.c means complex conjugate. To obtain the diffusion term, a formal but fast way to do consist in 
replacing the component of the test function Xj by ftp- in (^Kia^{w)'){X),X)^^ and (L(a^ (w)) (A), A)^^ 

defined by (99). In fact, in the proof of Proposition 6.7 one can see that the real part A] corresponds to 

i9yi and the imaginary A^ corresponds to d^2 , so that A, = X] + iX'j correspond to dyi + id^-^ = 2d^. 
3 3 j -J J J j j ± j 

That completes the proof of Proposition |6.3| ■ 
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